'  RD-A140  001 
UNCLASSIFIED 


DEVELOPMENT  OF  A  THIRD- INVARIANT  PLASTICITV  THEORV  FOR 
CONCRETE  AND  SOILS<U>  NEH  MEXICO  ENGINEERING  RESEARCH 
INST  ALBUQUERQUE  H  L  SCHREVER  FEB  84  AFWL-TR-82-119 
F29601-81-C-0012  F/G  20/11 


1/1 


NL 


AFWL-TR-83-1 19 


AFWL-TR- 

83-119 


DEVELOPMENT  OF  A  THIRD-INVARIANT  PLASTICITY 
THEORY  FOR  CONCRETE  AND  SOILS 

Howard  L.  Schreyer 

University  of  New  Mexico 

New  Mexico  Engineering  Research  Institute 

University  Station 

Albuquerque,  NM  87131 


February  1984 


Final  Report 


Approved  for  public  release;  distribution  unlimited. 


AIR  FORCE  WEAPONS  LABORATORY 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base,  NM  87117 


s 


DTIC 

ELECTE 
APR  1  1 1984 


84  04  io  023 


AFWL-TR-83-119 


This  final  report  was  prepared  by  the  New  Mexico  Engineering  Research 
Institute,  Albuquerque,  New  Mexico,  under  Contract  F29601-81-C-0013,  Job  Order 
2307 Y201  with  the  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  New 
Mexico.  Timothy  J.  Ross  (NTES)  was  the  Laboratory  Project  Officer-in-Charge. 

When  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  in  connection  with  a  definitely  Government-related  procure¬ 
ment,  the  United  States  Government  incurs  no  responsibility  or  any  obligation 
whatsoever.  The  fact  that  the  Government  may  have  formulated  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be 
regarded  by  implication,  or  otherwise  in  any  manner  construed,  as  licensing 
the  holder,  or  any  other  person  or  corporation;  or  as  conveying  any  rights  or 
permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in  any 
way  be  related  thereto. 

This  report  has  been  authored  by  a  contractor  of  the  United  States 
Government.  Accordingly,  the  United  States  Government  retains  a  nonexclusive, 
royalty-free  license  to  publish  or  reproduce  the  material  contained  herein,  or 
allow  others  to  do  so,  for  the  United  States  Government  purposes. 

If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing 
list,  or  if  your  organization  no  longer  employs  the  addressee,  please  notify 
AFWL/NTES,  Kirtland  AFB,  IW  87117  to  help  us  maintain  a  current  mailing  list. 

This  technical  report  has  been  reviewed  and  is  approved  for. publ i cation. 


PAUL  E.  MINTO 
Captain,  USAF 
Chief,  Applications  8ranch 


FOR  THE  COMMANDER 


STORM 
Monel ,  USAF 

lief.  Civil  Engineering  Research  Div 


00  NOT  RETURN  COPIES  OF  THIS  REPORT  UNLESS  CONTRACTUAL  OBLIGATIONS  OR  NOTICE 
ON  A  SPECIFIC  DOCUMENT  REQUIRES  THAT  IT  BE  RETURNED. 


REPORT  DOCUMENTATION  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 

1  ass  1  fled  _ 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  oeclassification/oownqraoing  schedule 


1b.  restrictive  markings 


3.  DISTRIBUTION/AVAILABILITY  of  report 

Approved  for  public  release;  distribution 
unlimited. 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER  IS) 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

AFWL-TR-83-11 9 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

University  of  New  Mexico 
New  Mexico  Engrg  Rsch  Inst 


6a  AOORESS  (City.  State  and  ZIP  Coda) 

University  Station 
Albuquerque,  NM  87131 


6a  NAME  OF  PUNOING/SPONSORING 
ORGANIZATION 


Sc.  AOORESS  (City,  State  and  ZIP  Coda) 


b.  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 
(IfappU cable) 

Air  Force  Weapons  Laboratory 


7b.  AOORESS  (City.  State  and  ZIP  Coda) 

Kirtland  Air  Force  Base,  NM  87117 


11.  TITLE  (Include  Security  damiftemtion!  DEVELOPMENT  OF  A 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F29601-81-C-001 3 


io.  source  of  funding  nos. 


PROGRAM 
ELEMENT  NO. 

611 02 F 


PROJECT 

TASK 

NO. 

NO. 

2307 

Y2 

WORK  UNIT 
NO. 


12.  PERSONAL  AUTHORtS) 

Schreyer,  Howard  L. 


13a  TYPE  OF  REPORT 

Final 


16.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED 
FROM  £2.  .QCt_  TO 


14.  OATE  OF  REPORT  (Yr..  Ho..  Dan) 

1984,  Februar 


IB.  PAGE  COUNT 

92 


COSAT  I  COOES 


GROUP 


11 


la  SUBJECT  TERMS  (Contlnua  on  tent  no  it  naeaaaary  and  identify  by  Mock  number) 

Viscoplasticity  Strain  Rate 

Concrete  Cracking 

Soils 


IB.  ABSTRACT  1  Continue  on  tea* tea  if  naeaaaary  and  identify  by  Mocb  number) 

A  viscoplastic  constitutive  relation  for  frictional  materials  has  been  formulated  and 
applied  to  concrete  and  soils.  First  and  third  invariants  of  both  stress  and  strain  are 
used  Instead  of  the  more  conventional  second  invariants  employed  in  relations  for  metal 
plasticity.  Detailed  comparisons  between  theoretical  and  experimental  data  are  made  for  a 
variety  of  paths  Involving  weak  concrete  and  sandy  soils.  For  the  soils,  both  static  and 
dynamic  behavior  are  considered.  The  model  has  been  extended  to  Incorporate  cracking  based 
on  a  maximum  principal  stress  criterion.  Proposed  experimental  and  theoretical  work  will 
cover  rate  effects  In  concrete,  the  response  of  unsaturated  clay,  and  Interfaces  between 
soil  and  concrete.. 


20.  OISTRISUTION/AVAI LABILITY  OF  ABSTRACT 

UNCLASSIP I EO/UN LIMITED  Q)  SAME  AS  RPT.  □  OTIC  USERS  □ 


22b.  NAME  OF  RESPONSIBLE  INOIVIOUAL 

Dr  Timothy  J.  Ross 


31.  ABSTRACT  SECURITY  CLASSIFICATION 


Unclassified 


22b.  TELEPHONE  NUMBER 
,  (Include  Area  Coda) 

(505)  844-9087 


22c.  OFFICE  SYMBOL 

NTES 


CONTENTS 


Section  Page 

I  INTRODUCTION  3 

II  RATE  EFFECTS  IN  CONCRETE  5 

Existing  Data  5 

Experimental  Device  for  Determining  Rate 

Effects  In  Concrete  11 

III  VISCOPLASTIC  CRACKING  MODEL  18 

IV  INTERFACE  MODELING  30 

V  RECOMMENDATIONS  32 

REFERENCES  34 

APPENDIXES 

A.  A  THIRD-INVARIANT  PLASTICITY  THEORY  FOR 

LOW-STRENGTH  CONCRETE  37 

8.  A  THIRD- INVARIANT  VISCOPLASTICITY  THEORY  FOR 

RATE-DEPENDENT  SOILS  65 


Accession  For 

NTIS  GRA&I 
DTIC  TAB 
Unannounced 
Just if Icat ton 


By - 

Distribution/ 


Availability  Codes 
Avail  and/or 
Dlst  |  Special 


1 


»  ■  ■  ilf. t J r.T'. y.V. 7.  V.V  A-.1  ;r*?V ”>  ■V/.-yA\T,MS.TC»7A^V.M 


ILLUSTRATIONS 

P23i 


Pressure  histories  for  some  propellants  14 

Conceptual  test  apparatus  for  tests  at  high 

strain  rates  15 

Behavior  of  cracking  function,  g,  with  application 
to  uniaxial  stress  20 

Plasticity  and  cracking  flow  surfaces  21 

Flow  chart  for  viscoplastic-cracking  algorithm  25 

Uniaxial  stress  26 

Cracking  failure  for  shear  stress  versus  shear  strain  27 

Indirect  shear  response  29 


y.% 


y 


v: 


I.  INTRODUCTION 


The  ability  to  predict  the  response  of  structures  subjected  to  large, 
abrupt  bursts  of  energy  depends  on  a  number  of  factors.  These  include  load 
definition,  the  characteristics  of  the  geological  material  In  which  the  struc¬ 
ture  Is  embedded,  the  mechanism  by  which  forces  are  transferred  from  the  geo¬ 
logical  material  to  the  structure,  and  the  characteristics  of  the  structure. 
The  basic  objective  of  the  research  project  discussed  in  this  report  is  to 
provide  an  improved  understanding  of  the  response  of  common  geological  and 
structural  materials  and  of  structure-media  (concrete-soil)  Interfaces  to 
abrupt  loading.  For  this  purpose,  engineering  constitutive  models  have  been 
developed,  and  appropriate  experimental  data  are  being  used  to  validate  the 
models. 

« 

The  Initial  formulation  of  a  constitutive  model  for  frictional  materials 
Involved  the  use  of  the  theory  of  plasticity  modified  to  Incorporate  first  and 
third  Invariants  of  stress  and  strain  Instead  of  the  more  conventional  second 
Invariants.  A  nonassoclated  flow  rule  was  applied  to  control  dilatation. 

This  theory  Is  outlined  In  Appendix  A,  where  detailed  comparisons  of  the  theo¬ 
retical  data  with  stress-strain  data  obtained  In  experiments  on  weak  concrete 
may  also  be  found.  The  experimental  stress  paths  encompass  a  wide  range  and 
provide  a  comprehensive  test  for  the  model. 

Rate  effects  represent  an  Important  feature  that  must  be  Included  If 
deformations  resulting  from  abruptly  applied  loads  are  to  be  predicted  accu¬ 
rately.  A  viscoplastic  formulation,  applied  to  sand  and  clayey  sands.  Is 
given  In  Appendix  B.  Theoretical  and  experimental  data  are  compared  for  both 
static  and  dynamic  loads,  and  It  Is  concluded  that  the  model  Is  a  reasonable 
base  for  representing  these  soil  materials  as  well  as  concrete. 

There  Is  a  dearth  of  experimental  data  on  multi  axial  rate  effects  for 
concrete  and  geological  materials.  A  survey  of  the  literature  (Section  II) 
has  disclosed  that  most  data  for  concrete  are  based  on  uniaxial  stress  condi¬ 
tions  with  a  limited  range  of  strain  rates.  Also  In  Section  II  Is  a  descrip¬ 
tion  of  a  conceptual  experimental  apparatus  In  which  propellants  will  be  used 
to  provide  forces  with  rise-times  sufficiently  small  that  multlaxlal  data  for 
concrete  with  strain  rates  on  the  order  of  1/s  can  be  obtained. 


The  viscoplastic  model  (Appendix  8)  has  been  modified  to  include  crack¬ 
ing.  The  procedure,  described  in  Section  III,  consists  of  using  a  maximum 
normal  tensile  stress  criterion  combined  with  strain  control  of  crack  propaga 
tion.  This  method  provides  nunerical  stability  together  with  a  physically 
realistic  model  that  can  also  be  used  with  an  interface  model. 

Several  interesting  developments  with  regard  to  the  representation  of 
interface  phenomena  have  been  reported  in  the  literature.  Some  of  the  impli¬ 
cations  of  these  studies  for  the  development  of  an  interface  model  for  soil 
and  concrete  are  described  in  Section  IV. 


The  primary  focus  for  the  continuation  of  this  research  project  is  summa 
rized  briefly  in  Section  V.  The  effort  will  Include  the  development  of  an 
experimental  device  for  obtaining  strain  rate  data  on  concrete,  the  develop¬ 
ment  of  interface  modeling,  and  both  theoretical  and  experimental  investiga¬ 
tions  of  the  response  of  wet  clay  to  tri axial  loading. 


II.  RATE  EFFECTS  IN  CONCRETE* 


m 


EXISTING  OATA 

Strain  rates  affect  several  of  the  mechanical  properties  of  concrete. 

The  extent  to  which  these  properties  are  affected  is  not  easily  determined 
because  it  depends  upon  the  property  being  considered  and  the  method  used  to 
prepare  and  test  the  concrete.  In  most  of  the  literature  It  Is  reported  that 
Increasing  the  strain  rate  at  which  a  test  Is  conducted  Increases  the  compres¬ 
sive  strength,  the  modulus  of  elasticity,  and  the  energy  absorption  of  the 
concrete.  The  sensitivity  of  a  concrete  to  dynamic  loads  is  related  to  the 
extent  of  the  Increase  of  these  parameters 

Testing  concrete  at  high  strain  rates  began  In  the  early  1900s  when  the 
effect  of  Impact  driving  on  concrete  piles  became  a  matter  of  Interest.  Since 
that  time  Interest  In  these  effects  has  grown,  and  the  tests  have  become 
Increasingly  sophisticated.  A  recent  literature  survey  indicates  this  growth 
(Refs.  1-6),  but  large  gaps  In  knowledge  are  also  evident.  Mainstone 
(Ref.  1)  tried  to  bring  together  the  available  data  on  the  sensitivity  to 
strain  rates  of  the  stiffness  and  nominal  compression  strength  of  concrete. 
However,  as  he  noted,  his  plot  of  the  ratio  of  dynamic  unconfined  compression 
strength,  f^,  to  static  unconfined  compression  strength,  f£,  versus  strain 

rate,  e,  does  not  take  Into  consideration  such  factors  as  aggregate  stiffness 
and  moisture  content.  Suarls  and  Shah  (Ref.  5),  on  the  other  hand,  have  tabu¬ 
lated  the  results  of  many  strain  rate  tests.  Including  tension,  flexure,  and 
compression  tests.  These  authors  note  that  concrete  Is  most  sensitive  to 
strain  rate  effects  In  tension  and  least  sensitive  In  compression.  They  do 
not  attempt  to  sort  out  the  parameters  that  affect  the  sensitivity  in  any  one 
mode  of  testing  (l.e.,  tension,  flexure,  compression). 

In  the  work  discussed  In  Refs.  1,  4,  5,  and  7,  the  first  result  obtained 
and  verified  was  that  as  6  Increases,  f^  also  Increases.  An  Immediate  second 
result  observed  was  that  the  amount  of  plastic  strain  decreases  as  strain 
rates  Increase  (Ref.  7).  In  1936,  Jones  and  Rlchart  (Ref.  8)  performed  some 
tests  on  concrete  aged  7  days  and  28  days.  The  fastest  strain  rate  at  which 

*Th1s  section  prepared  by  Susan  M.  Babcock,  Research  Engineer,  New  Mexico 
Engineering  Research  Institute,  University  of  New  Mexico,  Albuquerque. 


twy  tested  was  on  the  order  of  10-3/s.  Their  results  also  showed  that  an 
Increase  In  the  rate  of  loading  resulted  In  an  Increase  In  strength. 

Watsteln  (Ref.  9)  conducted  a  series  of  tests  on  weak  concrete  at  strain 
rates  of  about  0.004/s,  1.5/s  and  10/s,  and  on  higher  strength  concrete  at 
strain  rates  of  about  0.003/s,  2.9/s,  and  6.7/s.  An  hydraulic  machine  was 
used  in  performing  the  slower  tests  and  a  drophammer  in  the  faster  tests. 
Watsteln  concluded  that  the  compressive  strength,  the  modulus  of  elasticity, 
the  strain  energy  absorption,  and  the  strain  at  peak  stress  all  Increase  with 
Increasing  strain  rate. 

In  1958,  Evans  (Ref.  10)  reported  on  a  series  of  tests  performed  in  the 
1930s  In  which  a  compressed-air  testing  device  was  used.  He  concluded  that  no 
definite  Increase  in  strength  was  apparent  at  loading  speeds  to  failure  slower 
than  1/20  s.  At  faster  speeds  he  observed  a  definite  increase  In  the  crushing 
strength.  Evans  also  observed  that  the  percent  of  increase  In  crushing 
strength  with  Increase  in  load  rate  was  greater  for  the  leaner  mix  (high 
water -to -cement  ratio)  than  for  the  richer  mix  (low  water-to-cement  ratio)  of 
concrete. 

During  1964,  Green  (Ref.  3)  performed  some  tests  designed  to  determine  a 
specimen's  ability  to  absorb  energy  by  withstanding  repeated  blows  from  an 
impacting  pendulum  (Impact  strength).  Green  concluded  that  for  his  set  of 
specimens  (1)  the  Impact  strength  of  Portland  cement  concrete  increases  with 
the  compressive  strength;  (2)  the  impact  strength  Is  greater  for  concretes  In 
which  an  angular  rough  aggregate  Is  used;  (3)  the  impact  strength  Is  slightly 
less  for  water -cured  concrete  than  for  air-cured  concrete,  and  water-cured 
concrete  resists  cracking  better  but  fails  more  quickly  after  the  initial 

crack;  (4)  the  effects  of  sand  grading  are  slight;  and  (5)  the  modulus  of 
rupture  Is  not  a  good  indicator  of  impact  strength. 

In  1966  Cowell  (Ref.  11)  performed  a  series  of  tests  In  which  he  investi¬ 
gated  the  effects  of  compressive  strength  of  concrete  and  curing  technique  on 
strain  rate  effects.  He  obtained  strain  rates  of  from  about  10-6/s  to  0.8/s. 
Cowell  describes  thoroughly  his  testing  technique,  specimen  preparation,  and 
results.  He  used  medlun-  and  high-strength  concretes  and  cured  specimens  of 
both  strengths  In  two  ways  (1)  26-day  wet  cure  followed  by  2  days  at  room 


hunldlty  and  (2)  28-day  wet  cure  followed  by  21  days  at  room  humidity.  He 
observed  that  drying  the  concrete  tended  to  lower  the  rate  of  Increase  In 
compressive  strength  at  lower  stress  rates  for  both  strengths  or  concrete. 
However,  the  effects  of  drying  on  the  dynamic  strength  of  the  concrete  were 
more  pronounced  in  the  higher  strength  concrete.  He  also  noted  that  drying 
the  concrete  reduced  the  rate  of  Increase  of  the  secant  modulus  at  lower 
stress  rates  but  produced  negligible  effects  at  the  maximum  stress  rate. 

Cowell  also  observed  that  the  modulus  of  elasticity,  the  f'  value,  and 
Poisson's  ratio  ( v)  increased  with  increasing  loading  rates  for  the  stress 
rates  he  used.  When  he  computed  the  ratio  of  f^  tof^,  the  difference  between 
the  high-strength  and  medium-strength  concretes  was  Insignificant  and  could  be 
estimated  by  a  single  curve  In  the  space  of  f£/f£  versus  strain  rate. 

Also  in  1966,  Goldsmith,  Polivka,  and  Yang  (Ref.  12)  published  a  paper  in 
Experimental  Mechanics  describing  an  experimental  program  in  which  they  had 
used  the  Hopklnson  split  bar  technique.  Discussions  of  this  technique,  which 
was  originally  used  for  metals,  can  be  found  In  References  6,  13,  and  14.  The 
stress  levels  achieved  by  Goldsmith  et  al.  (Ref.  12)  were  too  low  to  produce 
significant  damage.  By  comparing  the  modulus  of  elasticity  from  shocked  por¬ 
tions  of  the  long  cylindrical  specimen  to  the  modulus  from  the  unshocked  por¬ 
tions,  they  were  able  to  conclude  that  the  dynamic  modulus  is  significantly 
greater  than  the  corresponding  static  modulus.  The  authors  also  studied  the 
dispersion  and  attenuation  of  the  waves  propagating  through  the  specimen. 

From  this  study  they  ascertained  that  the  longitudinal  waves  propagated  with¬ 
out  dispersion  and  with  little  attenuation. 

In  1967  Atchley  and  Furr  (Ref.  7)  performed  some  very  slow  (virtually 
static)  tests  and  some  drophammer  tests  on  concrete.  The  authors  concluded 
that  concrete  can  absorb  only  a  limited  amount  of  energy,  after  which  point 
failure  will  occur.  However,  they  found  It  difficult  to  determine  the  cause 
of  the  initial  failure.  They  computed  values  of  fyf£  for  their  data  and 
compared  these  results  with  data  from  Cowell  (Ref.  11)  and  Watsteln  (Ref.  9). 
In  the  space  of  f^/f£  versus  A,  both  Cowell's  and  Watsteln' s  data  formed  con¬ 
cave  upward  curves;  Atchley's  and  Furr's  data  formed  concave  downward  curves. 
Atchley  and  Furr  tentatively  attributed  the  discrepancy  to  a  difference  In 


strain  gage  locations.  Cowell  and  Watsteln  had  both  used  surface  strain 
gages;  Atchley  and  Furr  had  used  Internal  strain  gages.  However,  other  dif¬ 
ferences  among  these  three  research  efforts  may  have  also  contributed  to  the 
differing  characteristics  of  the  curves.  Among  these  factors  were  differences 
in  curing  techniques,  concrete  age  at  test  time,  and  type  of  coarse  aggregate 
used.  Atchley  and  Furr  did  concur  with  previous  researchers  that  concrete 
strength  Increases  with  Increasing  strain  rate  but  seems  to  approach  asynptot- 
ically  some  undetermined  limit  within  the  range  of  strain  rates  measured. 

In  the  1970s  and  early  1980s,  several  papers  were  published  in  which 
parametric  studies  at  high  strain  rates  were  discussed  (Refs.  2,  4,  and  15- 
19).  The  parameters  studied  included  cement/water  ratio  (Refs.  8  and  16),  age 
(Refs.  2  and  16),  cement  content  (Refs.  8  and  16),  aggregate  stiffness 
(Refs.  16-18),  coarse/fine  aggregate  ratio  (Ref.  16),  curing  conditions 
(Refs.  2,  4,  and  15),  and  compressive  strength  (Ref.  4).  According  to  Hughes 
and  Gregory  (Ref.  16)  and  Hughes  and  Watson  (Ref.  18),  the  dynamic  strength  of 
concrete  Increases  slightly  more  at  a  high  water/cement  ratio  than  at  lower 
ratios.  The  effect  of  the  cement  volume  fraction  on  the  dynamic-strength/ 
static-strength  ratio  appears  to  have  been  insignificant  in  these  two  studies, 
and  the  authors  also  report  that  variations  in  type  of  coarse  aggregate  used 
had  Insignificant  effects  on  the  dynamic  strength  of  the  concrete.  In  one 
case  the  aggregates  were  limestone  and  gravel  (Ref.  8)  and  in  the  other  case, 
limestone  and  granite  (Ref.  16).  Sparks  and  Menzies  (Ref.  17),  on  the  other 
hand,  concluded  that  the  "sensitivity  of  the  static  strength  to  the  rate  of 
loading  appears  to  be  related  to  the  stiffness  of  the  aggregate."  Their  work 
showed  that  stiff  aggregates  such  as  limestone  provided  less  Increase  in 
strength  with  Increased  load  rate  than  did  less  stiff  aggregates  such  as  gran¬ 
ite.  The  least  stiff  aggregate  they  tested  gave  the  greatest  Increase  in 
strength  with  increasing  load  rate. 

The  coarse-to-f ine  aggregate  ratio  results  are  plotted  in  Reference  16. 
The  plot  Indicates  that  the  reduction  in  strength  effected  by  the  lower 
water/cement  ratio  was  not  as  great  when  a  larger  amount  of  fine  aggregate  was 
used.  The  effects  of  the  age  of  the  concrete  on  the  dynamic  strength,  as 
reported  by  Hughes  and  Gregory  (Ref.  16),  were  Inconsistent.  Conversely, 
Kaplan  (Ref.  2)  observed  that  as  the  age  of  the  concrete  increased  from  2  days 


to  29  days,  the  rate  of  increase  of  the  dynamic  strength  with  increasing  load 
rate  was  more  linear  at  load  rates  from  0.001  to  100  N/mm2s.  The  data  for  the 
newer  concrete  tended  to  have  a  definite  concave  upward  curve  in  the  space  of 
compressive  strength  versus  rate  of  loading. 

The  curing  condition  used  for  the  concrete  has  been  observed  to  have  a 
significant  effect  on  the  sensitivity  of  the  concrete  to  high  strain  rates 
(Refs.  2  and  4).  Ohir  and  Sangha  (Ref.  4)  compared  data  found  in  the  litera¬ 
ture  (Refs.  8,  9,  and  10)  as  well  as  data  they  had  obtained.  They  concluded 
that  as  the  length  of  the  air-cure  period  is  increased,  the  effect  of  strain 
rate  diminishes.  This  finding  is  supported  by  Kaplan's  observations  (Ref.  2) 
that  moist  specimens  were  significantly  more  sensitive  to  increasing  strain 
rates  than  were  drier  specimens.  Spooner  (Ref.  15)  tested  for  sensitivity  of 
strain  rate  to  cure  method  for  relatively  low  strain  rates  (3.8  x  10_7/s  and 
5.0  x  10”5/s).  His  finding  that  air  cure  decreases  the  limit  stress  more  than 
does  water  cure  agrees  with  the  results  presented  in  References  2  and  4.  He 
also  noted  that  strain  at  the  maximum  stress  does  not  differ  significantly  on 
the  two  curves.  Kaplan  (Ref.  2)  suggests  that  the  enhanced  dynamic  strength 
of  moist  concrete,  as  compared  to  dry  concrete,  is  due  to  the  hydrostatic 
pressure  in  the  pores,  which  delays  microcracking  and  allows  a  higher  compres¬ 
sive  strength  to  develop. 

Ohir  and  Sangha  (Ref.  4)  also  studied  the  effect  of  the  nominal  strength 
of  the  concrete  on  the  sensitivity  of  the  concrete  to  increasing  strain  rates, 
as  reported  in  the  literature.  They  found  conflicting  data.  They  reported 
that  some  data,  such  as  those  obtained  by  Watstein  (Ref.  9),  showed  a  slight 
Increase  In  sensitivity  to  higher  strain  rates  when  the  nominal  strength  of 
the  concrete  was  increased.  This  reported  trend  in  Watstein' s  data  does  not 
seem  to  be  confirmed  by  Watstein's  own  report.  Other  data  (Refs.  8  and  10) 
Indicated  that  sensitivity  to  strain  rate  Is  Independent  of  concrete  strength. 
On  the  basis  of  their  own  experience  with  sandstones  ("natural"  concretes)  at 
relatively  low  rates  of  loading,  Ohir  and  Shangha  (Ref.  4)  concluded  that  the 
percent  increase  in  the  dynamic  strength  of  higher  strength  concretes  is 
greater  than  that  of  low-strength  concretes.  Some  discrepancies  regarding 
the  enhancement  of  dynamic  concrete  strength  and  its  association  with  static 
concrete  strength  are  evident  in  the  findings  of  various  authors.  The 
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conclusions  of  Dhlr  and  Sangha,  for  example,  differ  from  the  conclusions  Main- 
stone  (Ref.  1)  reached  after  reviewing  several  papers  (among  others.  Refs.  8, 
9,  10,  16,  and  17). 

For  the  data  discussed  so  far,  the  compression  tests  were  uniaxial  stress 
tests.  Only  two  types  of  multlaxial  tests  were  described  in  the  literature. 
Biaxial  testing  with  a  very  large  hydropneumatic  biaxial  device  and  a  device 
based  on  the  Hopkinson  bar  has  been  accomplished  on  metals  and  has  been  pro¬ 
posed  for  concrete  (Ref.  19).  The  other  type  of  multiaxial  test  was  performed 
by  Takeda,  Tachikawa,  and  Fijimoto  (Ref.  20).  These  workers  conducted  some 
trlaxial  tests  in  which  the  confining  pressure  was  applied  statically  and  the 
axial  stress  was  applied  at  a  high  strain  rate.  They  drew  no  conclusions 
regarding  multlaxial  rate  effects. 

On  the  basis  of  the  data  obtained,  the  following  conclusions  can  be 
drawn: 

1.  As  the  loading  or  strain  rate  increases,  fj  increases.  There  is  some 
dispute  about  whether  the  ratio  of  f^/f£  continues  in  an  accelerating  fashion 
when  plotted  against  e,  as  concluded  in  Refs.  4,  9,  10,  and  11  among  others, 
or  levels  off  as  e  becomes  very  large  (Ref.  7). 

2.  As  the  strain  rate  increases,  the  amount  of  plastic  strain  decreases. 
A  related  factor  is  the  increase  found  in  the  secant  modulus  and  the  late-time 
tangent  modulus.  The  initial  tangent  modulus  apparently  remains  constant  over 
the  spectrum  of  strain  rates. 

3.  The  type  of  coarse  aggregate  used  in  the  concrete  has  a  major  effect 
upon  the  sensitivity  of  the  concrete  to  strain  rate  effects.  The  stiffer 
aggregates  give  less  sensitivity  than  do  the  less  stiff  aggregates. 

4.  The  moisture  content  of  the  concrete  has  a  significant  influence  on 
the  sensitivity  of  the  concrete  to  strain  rate  effects.  Concretes  with  high 
moisture  contents  are  more  sensitive  than  are  drier  concretes. 

A  data  base  is  beginning  to  emerge  for  uniaxial  compression  tests  on 
concrete  at  a  variety  of  strain  rates.  Even  though  some  of  the  data  conflict, 
one-dimensional  constitutive  models  can  be  developed  by  attempting  to  match 
the  most  relevant  data  set  for  a  particular  application.  To  expand  to  more 


10 


dimensions,  assumptions  and  engineering  judgment  must  be  used  because  the 
are  virtually  no  data  for  multiaxlal  strain  rate  effects.  Bazant  and  Oh 
(Ref.  21)  commented  that  three-dimensional  effects  such  as  inelastic  compac¬ 
tion,  shear  enhanced  compaction,  and  dilatation  probably  are  not  significant 
at  high  strain  rates  and  can  be  neglected  in  the  development  of  models.  On 
the  other  hand,  these  three-dimensional  effects  have  only  recently  been  recog¬ 
nized  for  their  Importance  in  the  static  three-dimensional  modeling  of  con¬ 
crete.  Therefore,  to  determine  the  Importance  of  three-dimensional  effects  at 
high  strain  rates,  multiaxial  high-strain-rate  tests  must  be  performed.  From 
the  uniaxial  stress  tests,  significant  parameters  are  becoming  apparent. 

These  same  parameters  may  prove  to  be  Important  in  the  multiaxial  tests  dis¬ 
cussed  in  the  next  subsection.  On  the  basis  of  the  papers  reviewed.  It  is 
clear  that  the  parameters  and  techniques  used  in  the  proposed  multiaxial  tests 
must  be  well  defined  and  documented.  This  matter  will  become  particularly 
Important  as  later  researchers  try  to  Incorporate  the  test  data  Into  a  unified 
data  base. 

EXPERIMENTAL  DEVICE  FOR  DETERMINING  RATE  EFFECTS  IN  CONCRETE 

Data  from  multiaxial  stress  tests  at  high  rates  of  strain  are  virtually 
nonexistent.  Because  of  the  lack  of  multiaxial  data,  any  attempt  to  model  the 
multiaxial  behavior  of  concrete  requires  that  assunptlons  be  made.  The 
effects  of  high  strain  rates  on  the  three-dimensional  behavior  of  concrete  are 
unknown.  To  fill  this  void,  a  relatively  simple  and  inexpensive  testing  pro¬ 
cedure  is  proposed. 

Concrete  is  used  in  a  number  of  applications  that  expose  the  concrete  to 
high  strain  rates.  Airblast  and  ground  motion  loadings  from  explosive  shocks 
are  two  examples  of  such  exposure.  Other  examples  include  the  impact  driving 
of  concrete  piles,  the  accidental  dropping  of  concrete  structural  members 
during  handling,  high  winds,  earthquakes,  and  sudden  increases  in  pressure  in 
concrete  vessels.  Most  of  these  events  can  Involve  strain  rates  ranging  from 
virtually  static  (10"6/s)  to  as  high  as  106/s.  It  has  been  shown  that  at 
strain  rates  of  less  than  about  0.1/s,  insignificant  or  confusing  and  uncon- 
flrmable  results  are  obtained  (Ref.  1).  On  the  other  hand,  for  the  large 
majority  of  strain  environments  that  protective  structures  are  expected  to 
survive,  particularly  airblast,  ground  shock,  wind,  and  earthquake,  the  strain 


rates  are  less  than  10/s.  For  these  reasons,  the  strain  rates  sought  in  an 
experimental  program  should  fall  between  0.1/s  and  10/s. 

The  rise-time  to  peak  stress  can  be  determined  for  a  particular  concrete 
once  e  is  given.  The  stress  rate,  a,  can  be  approximated  for  uniaxial  stress 
by 

i  *  Ee 

where  E  is  the  modulus  of  elasticity.  To  find  the  approximate  rise-time  for  a 
given  strain  rate,  the  stress  rate  is  calculated,  and  the  time  required  to 
reach  f '  can  then  be  determined.  For  example,  at  a  strain  rate  of  10/s  and  a 
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compressive  strength  of  27.6  MPa  (an  average-strength  concrete),  E  is  calcula¬ 
ted  to  be  2.5  x  104  MPa  (Ref.  22),  and  o  Is  25  x  104  MPa/s.  Therefore,  the 
rise-time  Is  26.7  MPa/25  x  104  MPa-s,  or  approximately  0.1  ms. 

The  size  of  the  concrete  test  specimen  is  constrained  by  two  factors. 

The  minimum  size  is  limited  by  the  nonhomogeneous  nature  of  concrete.  The 
specimen  must  be  large  enough  to  ensure  that  the  nonhomogeneities  in  the  con¬ 
crete  can  be  neglected.  To  obtain  this  homogenization,  it  is  generally 
assumed  that  the  specimen  must  be  at  least  10  times  the  size  of  the  largest 
aggregate  In  the  mix.  Usually,  a  3/8-in  coarse  aggregate  is  used,  which  means 
that  the  minimum  specimen  size  Is  about  100  mm.  Conversely,  to  eliminate  the 
effects  of  Inertia  associated  with  the  passage  of  waves,  the  maximum  specimen 
size  must  be  limited  in  order  to  minimize  the  transit  time  of  the  wave  (or  to 
maximize  the  number  of  wave  reflections).  If  the  number  of  wave  reflections 
is  not  sufficiently  large  during  loading,  the  strain  across  the  specimen  can¬ 
not  be  considered  uniform  and  the  equations  of  motion  will  have  to  be  solved 
and  constitutive  Information  backed  out  indirectly.  For  a  large  number  of 
reflections  the  response  can  be  considered  homogeneous  and  quasi-static.  The 
transit  time  Is  calculated  from  the  wave  speed  in  concrete  (about  3000  m/s) 
and  from  the  specimen  size.  The  minimum  specimen  size  of  100  mm  gives  a  tran¬ 
sit  time  of  about  33  us.  With  a  strain  rate  of  10/s,  the  rise-time  is  calcu¬ 
lated  to  be  0.1  ms,  and  three  reflections  will  occur  before  a  stress  of  f*  Is 
reached.  A  strain  rate  of  1.0/s  requires  a  rise-time  of  10  ms,  and  300 


reflections  will  occur  before  f'  Is  reached.  Three  hundred  reflections  are 
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more  than  enough  to  ensure  that  the  effects  of  Inertia  will  not  significantly 
affect  the  results;  three  reflections  probably  are  not  enough.  The  minimum 
number  of  reflections  will  have  to  be  determined  analytically  by  a  comparison 
of  the  dynamic  solution  with  the  static  solution.  As  an  initial  estimate,  at 
least  25  to  50  reflections  will  be  required,  which  would  limit  the  highest 
quasi-static  6  to  about  1/s.  If  the  specimen  is  loaded  from  both  sides,  which 
effectively  reduces  the  transit  time  by  one-half,  a  quasi-static  6  of  about 
2/s  could  be  obtained.  Higher  strain  rates  can  be  obtained,  but  the  calcula¬ 
tion  of  stress  may  Involve  the  use  of  the  equation  of  motion. 

Figure  1  shows  typical  examples  of  pressure  histories  for  some 
propellants.  It  Is  evident  that  rise-times  of  less  than  10  ms  can  be 
achieved;  thus,  these  materials  could  be  used  to  obtain  the  required  loading 
rates.  Gunpowder,  another  readily  available  propellant,  can  produce  a  faster 
rise-time  as  well  as  higher  pressures  than  those  shown  In  Figure  1  and  thus 
represents  an  option  to  be  considered. 

The  proposed  testing  apparatus,  shown  In  Figure  2a,  Is  basically  a  thick- 
walled  steel  tube  bent  Into  an  oval.  The  two  pressure  chambers  are  shown  In 
Figure  2b.  The  high-pressure  chamber  maintains  a  pressure  at  which  the  pro¬ 
pellant  will  burn  properly;  the  low-pressure  chamber  provides  the  forcing 
function  for  the  specimen.  Variations  In  the  orifice  configuration  of  the 
steel  plate  separating  the  high-  and  low-pressure  chambers  can  be  used  to 
tailor  the  waveform.  Larger  holes  in  the  plate  will  produce  a  faster  rise¬ 
time  In  the  low-pressure  chamber,  and  smaller  holes  will  produce  a  slower 
rise-time.  The  exact  configurations  that  produce  the  desired  strain  rates  In 
the  concrete  must  be  determined  by  experimentation.  When  strain  gages  are 
mounted  on  a  platen,  the  platen  becomes  a  load  cell  that  can  be  used  with 
strain  gages  1  c  .  ed  on  the  concrete  specimens  to  obtain  stress-strain 
curves. 

The  concre  c  >tcimen  can  be  cylindrical  for  the  Initial  feasibility 
experiments.  Uniaxial  stress  tests  are  sufficient  for  this  preliminary  phase. 
When  cylindrical  specimens  are  used,  the  design  of  the  steel  platen  Is  con¬ 
siderably  simplified  because  the  platens  can  be  simple  cylinders  also.  Once 


Figure  1.  Pressure  histories  for  some  propellants 


Conceptual  test  apparatus  for  tests  at  high  strain  rates 


:eel  containment  vessel 


Figure  2.  Concluded 


the  apparatus  and  the  technique  have  been  thoroughly  tested,  modified  as  nec¬ 
essary,  and  verified,  more  complex  stress  path  tests  and  specimen  shapes  can 
be  used.  Cubical  specimens  can  be  used  to  advantage  in  multi  axial  stress  path 
tests  because  they  allow  a  variety  of  stress  paths.  However,  a  cubical  speci¬ 
men  requires  a  more  complicated  platen  design:  the  face  of  the  platen  must  be 
square,  and  the  back  must  be  circular  to  enable  the  platen  to  fit  Inside  the 
pressure  vessel. 

With  the  proposed  testing  apparatus,  a  great  variety  of  stress-defined 
path  tests  could  be  conducted.  Conventional  static  uniaxial  stress  tests  can 
be  performed  If  a  powerful  enough  air  compressor  can  be  found  or  If  very  slow 
burning  propellants  are  used.  Trlaxlal  tests  on  a  cylindrical  specimen  can  be 
undertaken  by  pressurization  of  the  lateral  surface.  Biaxial  static  or  quasi¬ 
static  tests  can  be  performed  on  cubical  specimens  If  two  sets  of  apparatuses 
on  orthogonal  axes  are  used.  In  any  of  these  modes,  strain  rate  tests,  at 
various  strain  rates,  can  be  accomplished  by  using  the  appropriate  amounts  of 
propellant. 


III.  VISCOPLASTIC  CRACKING  MODEL 


An  Important  part  of  the  research  effort  has  been  the  formulation  and 
development  of  an  algorithm  that  combines  cracking  with  the  conventional  plas¬ 
ticity  algorithm  for  frictional  materials.  Almost  any  problem  of  engineering 
significance  Involves  zones  In  which  one  or  more  principal  components  of  the 
stress  tensor  go  into  tension.  To  Ignore  this  feature  or  to  use  a  simple 
mean-pressure  cutoff  provides  misleading  results.  On  the  other  hand,  a 
detailed  cracking  algorithm  Is  often  more  accurate  than  warranted  by  the  needs 
of  the  analyst,  and  the  result  Is  an  Inefficient  numerical  procedure.  The 
approach  described  here  is  an  attempt  to  provide  an  approximation  procedure 
that  Is  both  fast  and  reasonably  accurate. 

Tensile  stresses  appear  In  a  variety  of  engineering  problems.  For  those 
Involving  frictional  materials  rather  than  metals,  the  situation  Is  more  com¬ 
plicated  because  frictional  materials  behave  differently  In  tension  and  In 
compression.  Consequently,  both  the  theory  and  the  computational  algorithm 
must  be  adjusted  to  reflect  actual  tensile  behavior.  This  behavior  Is  par¬ 
tially  understood  for  uniaxial  tensile  stress  only;  for  more  complex  states  of 
stress,  considerable  engineering  judgment  must  be  used.  Because  of  this 
uncertainty,  theoretical  and  numerical  algorithms  should  be  kept  as  simple  as 
possible  so  that  the  highest  possible  efficiency  and  also  a  level  of  accuracy 
consistent  with  the  available  experimental  data  can  be  achieved.  In  the  pro¬ 
cedure  proposed  here,  the  plasticity  algorithm  used  for  conventional  plastic¬ 
ity  In  the  compressive  regime  Is  used  to  simulate  cracking. 

Tensile  cut-off  behavior  Is  often  based  on  mean  pressure,  but  numerous 
results  Indicate  that  this  approximation  Is  not  adequate.  The  reason  Is  that 
one  stress  component  can  be  large  In  tension,  but  as  long  as  the  other  princi¬ 
pal  stresses  are  sufficiently  large  In  compression,  such  a  model  does  not 
predict  cracking.  Here,  It  Is  assumed  that  physical  behavior  Is  better  repre¬ 
sented  by  the  monitoring  of  Individual  components. 

Suppose  one  component  of  the  stress  tensor  goes  into  tension.  It  Is 
assumed  that  the  response  Is  elastic  up  to  the  cracking  stress,  o ,  and  that 
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cracking  does  not  occur  abruptly.  Instead,  cracking  occurs  gradually  as  a 
function  of  tensile  strain.  The  difference  between  total  strain  and  elastic 
strain  Is  called  the  cracking  strain  (ec).  The  cracking  stress  Is  also 
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assimed  to  reduce  linearly  with  total  tensile  strain.  The  rate  of  reduction 
is  large  but  finite,  and  this  finiteness  helps  to  stabilize  the  nunerical 
algorithm.  Also,  the  cracking  strain  limit,  e£,  is  a  parameter  that  can  be 
used  to  relate  energy  dissipation  to  finite-element  size  so  that  results  for 
crack  propagation  can  be  made  invariant  with  respect  to  the  finite-element 
mesh. 

To  accommodate  this  feature  within  a  plasticity  algorithm,  a  cracking 
potential  function  is  defined: 

♦c  ■  o  -  9(eC)  (1) 

where  o  (positive  In  tension)  denotes  the  stress  component  of  interest,  g  the 
cracking  function,  and  e  the  corresponding  cracking  strain  component.  Simi¬ 
lar  to  the  sign  convention  for  $  used  with  plasticity,  the  cracking  algorithm 
is  not  activated  if  $c  <  0  or  o  <  g.  Also,  $c  >  0  is  not  allowed,  and  crack¬ 
ing  Is  considered  to  be  occurring  if  ♦  *  0  or  c  «  g.  The  cracking  function, 

c  c  r 

g,  reduces  linearly  from  o0  to  zero  as  the  cracking  strain  Increases  to  e£, 

which  denotes  a  complete  fracture.  For  ec  >  e£,  g  is  zero. 

Also  in  analogy  to  conventional  plasticity,  the  total  tensile  strain  Is 
taken  to  be  the  sum  of  the  cracking  strain  and  the  elastic  strain: 

e  *  ec  +  ee  (2) 

and  the  stress  is  related  to  the  elastic  strain  through  the  conventional  elas¬ 
ticity  constitutive  relations.  If  a  stress  reversal  occurs,  the  crack  closes 
and  the  conventional  relations  for  plasticity  are  used  In  compression.  How¬ 
ever,  a  second  reversal  results  In  a  limited  tensile  stress  component  based  on 
the  accumulated  cracking  strain  (Fig.  3). 

To  Indicate  how  the  procedure  works  in  higher  dimensions,  consider  a  case 
of  plane  stress  (principal  stresses)  with  flow  surfaces  as  shown  in  Figure  4. 
The  solid  line  labeled  ♦  *  0  represents  a  conventional  plasticity  flow  sur¬ 
face,  while  and  $2  “  0  represent  cracking  flow  surfaces  for  stress 


19 


Figure  4.  Plasticity  and  cracking  flow  surfaces. 

components  and  o2,  respectively.  If  de{  and  de2  (positive  in  tension) 
denote  normal  components  of  inelastic  strain,  then  corresponding  components  of 
cracking  strain  are  defined  as  follows: 

de*  »  def  if  dej  >0  de£  *  de2  if  de2  >  0 


if  dej  <  0 


-  0  if  del  <  0 


«i  -  f<*e<i 


e2  *  fde\ 


That  is,  the  cracking  strains  are  considered  to  be  the  monotonic  accumulation 
of  Inelastic  tensile  Increments.  The  cracking  strains  that  are  obtained 
directly  from  surfaces  <^1  *  0  and  $2  *  0  are  not  assumed  to  affect  the 
plasticity  region.  However,  positive  Inelastic  strains  (dilatation)  from  the 
plasticity  surface  are  assumed  to  affect  the  cracking  surface.  A  detailed 
Investigation  of  experimental  data  will  be  required  to  determine  whether  these 
assumptions  are  reasonable. 


On  the  flow  surfaces  $!  *  0  and  $2  *  0,  an  associated  flow  rule  is  used. 
For  example,  consider  the  surface  limiting  the  first  component  of  stress: 

$1  *  -  9 i  (e?)  (5) 

where  g!  is  a  function  similar  to  the  function  g  shown  in  Figure  3.  The  ine¬ 
lastic  strain  increments  are  then  given  by  the  conventional  flow  rule  of 
plasticity: 


del  =  dX  iii 


(6) 


de2  a  dx 


111.  0 

3o2 


The  magnitude  of  the  scalar  parameter,  dx,  is  determined  by  enforcing  the  con¬ 
sistency  condition  3  0.  Thus,  for  excursions  in  stress  dominated  by  a 
tensile  component,  olt  the  inelastic  strain  component  dej  accumulates,  but  e2 
does  not  change  because  del  3  0*  As  cracking  continues,  ox  decreases  to  zero, 
representing  a  complete  crack,  and  «2  adjusts  automatically  according  to  the 
stress-strain  relations  of  elasticity. 

It  Is  assumed  that  the  cracking  strength  in  the  second  coordinate  direc¬ 
tion  is  not  affected  by  cracking  in  the  first  coordinate  direction  (cracks 
perpendicular  to  the  first  coordinate  direction)  and  the  reverse.  This 
assumption  is  carried  even  farther  for  the  case  of  shear  where  shear  failure 
is  assumed  to  be  Independent  of  tensile  cracking,  i.e.,  there  is  no  interac¬ 
tion  between  normal  and  shear  strains.  Although  it  is  known  that  the  assump¬ 
tion  cannot  be  completely  valid.  It  should  be  reasonable  for  materials  such  as 
concrete  because  of  aggregate  Interlock.  There  is  some  interaction,  but  only 
If  the  tensile  strains  become  large.  The  cracking  function  $12  used  for  shear 
is  similar  to  that  used  for  a  normal  component;  thus,  the  corresponding  crack¬ 
ing  function  (potential  function)  and  flow  relation  are 

$12  *  I012I  "9l2(le12l )  (8) 
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■  dx  sign  (<j12)  (9) 

where  g12  is  also  similar  to  g. 

The  problem  with  this  formulation  is  that  cracking  strains  can  accumulate 
for  states  with  a  large  mean  pressure,  which  is  not  realistic.  A  way  to  cir¬ 
cumvent  this  problem  is  to  scale  the  cracking  strain  increments  with  mean 
pressure,  just  as  the  inelastic  strain  invariant  is  scaled  with  mean  pressure. 

At  the  end  of  a  problem,  the  accumulated  cracking  strain  components  are 
available  at  each  node  or  element  according  to  the  numerical  formulation. 

These  tensor  components  can  be  expressed  in  terms  of  principal  values,  which 
are  associated  with  principal  directions.  If  a  principal  value  is  positive, 
an  appropriate  Interpretation  would  be  that  a  crack  has  formed  with  a  normal 
in  the  corresponding  principal  direction.  For  example,  if  only  a  shear  crack¬ 
ing  strain  is  developed,  this  could  correspond  to  a  crack  oriented  at  45  deg 
to  the  given  coordinate  axes. 

For  the  numerical  algorithm,  strain  increments  are  assumed  given  and  a 
trial  state  of  stress  based  on  an  elastic  step  is  obtained.  The  cracking 
functions  are  checked  first,  and  if  any  of  them  has  a  value  greater  than  zero, 
the  material  is  allowed  to  crack  until  a  redistribution  of  stress  and  strain 
reduces  the  value  to  zero.  The  plasticity  function  is  then  checked;  if  it  has 
a  value  greater  than  zero,  plastic  deformation  is  allowed  to  occur  until  this 
value  is  reduced  to  zero.  The  process  is  repeated  iteratively  until  all  func¬ 
tions  are  <  0.  The  algorithm  is  then  exited. 

This  process  of  satisfying  two  or  more  flow  functions  simultaneously  is 
actually  an  application  of  the  condition  that  at  a  corner  of  a  yield  surface 
that  does  not  have  a  continuous  slope,  the  inelastic  strain  increment  must 
fall  between  the  normals  of  the  respective  yield  functions  at  the  Intersection 
point.  All  equations  are  supposed  to  be  satisfied  simultaneously,  but  given 
the  nature  of  the  assumptions  on  cracking,  the  iterative  procedure  should  not 
introduce  unacceptable  errors. 


A  flow  chart  of  the  algorithm  is  given  in  Figure  5.  The  flow  functions 
are  defined  as  follows  for  states  of  stress  as  general  as  plane  stress  and 
plane  strain: 

»  L  -  yP  -  o$  (plasticity) 


*2  *  ai  ‘  9 1  ( e l ) 


(cracking  1-direction) 


<f>3  3  °2  "  92(e2) 


(cracking  2-direction) 


♦4  *  °3  *  93(e3) 


(cracking  3-direction) 


♦5  *  lo12l  -  g12(e3)  (cracking  1-2  shear) 


In  which  the  terms  In  4^  are  defined  In  Appendix  A;  o2,  and  o3  are  normal 
components  of  stress;  a12  is  a  shear  component;  and  the  functions  gL,  g2,  93* 
and  g12  are  similar  to  g. 


A  few  elementary  path  tests  were  used  to  check  the  basic  alc»,"'i  chm  in 
order  to  ensure  that  essential  response  characteristics  were  being  reproduced. 
The  paths  and  results  are  described  as  follows. 


Path  1,  uniaxial  stress  (Figs.  6a  and  6b) 

a.  Load  ou  into  compression. 

b.  Release  and  load  into  tension. 

c.  Continue  to  Impose  tensile  strain. 

d.  Load  Into  compression. 

e.  Superimpose  Increment  of  tensile  strain  (material  can  no  longer  sus¬ 
tain  any  tensile  stress). 


Path  2,  pure  shear  (Fig.  7) 

A  shear  strain  Is  applied  monotonical 1y,  which  causes  the  shear  stress  to 
build  to  the  cracking  level,  decrease  to  zero,  and  remain  there. 


vl scopl asti c-cracki ng  al gori thm. 


Path  3,  Indirect  shear  (Figs.  8a  and  8b) 

Here  the  first  component  of  strain  Is  decreased  while  the  second  compo¬ 
nent  Is  Increased  by  the  same  amount.  The  tensile  component  of  stress  dis¬ 
plays  the  cracking  phenomenon,  which  causes  a  simultaneous  adjustment  In  the 
compressive  stress  when  the  stress  redistribution  occurs. 

Although  these  examples  represent  elementary  stress  paths,  they  show  that 
the  algorithm  Is  operational.  The  theory  and  numerical  procedure  can  now  be 
used  as  the  basis  for  an  In-depth  engineering  model  Investigation  of  the 
cracking  of  concrete.  This  investigation  will  provide  guidelines  to  be  used 
in  seeking  suitable  experimental  data. 

Another  important  application  involves  interface  models  in  which  the 
possibility  of  separation  is  generally  not  considered.  Instead,  such  models 
attempt  to  predict  Interface  shear  stress  as  a  function  of  normal  pressure  and 
relative  motion.  For  dynamic  situations,  separation  is  a  distinct  possibil¬ 
ity.  Therefore,  a  cracking  algorithm  must  be  combined  with  an  interface  model 
to  provide  a  complete  analytical  capability. 


IV.  INTERFACE  MODELING 


Coulomb's  law  describes  frictional  effects  between  rigid  bodies  subject 
to  the  gross  sliding  of  one  body  relative  to  another.  Because  of  its  apparent 
conceptual  simplicity,  many  engineers  have  attempted  to  apply  the  law  point- 
wise  to  contact  problems  in  contlnuun  mechanics.  Even  for  linearly  elastic 
bodies,  there  are  fundamental  questions  concerning  the  existence  of  solutions 
(Ref.  23).  This  lack  of  a  proof  of  existence  may  be  one  of  the  reasons  that 
no  single  frictional  law  has  been  generally  accepted  as  providing  results 
consistent  with  observations.  For  metallic  surfaces,  a  nonlocal  and  nonlinear 
friction  law  is  suggested  (Ref.  23).  It  is  possible  that  a  similar  law  may  be 
appropriate  for  interfaces  between  concrete  and  soils. 

The  study  of  interface  phenomena  is  fraught  with  difficulties  as  indi¬ 
cated  by  the  results  presented  in  Reference  24.  When  a  constant  force  is 
applied  to  one  body  in  a  direction  parallel  to  an  interface,  a  constant  slip 
rate  is  often  observed.  However,  in  other  cases  that  appear  similar,  a  non¬ 
constant  slip  motion  (stick-slip)  is  observed.  It  can  be  shown  that  this 
instability  can  occur  if  the  stiffness  of  the  moving  body  is  less  than  a  cri¬ 
tical  value  that  depends  on  the  inertia  of  the  body. 

The  possibilities  of  nonexistent  solutions  and  instabilities  are  rarely 
considered  in  engineering  models  of  interfaces.  Some  models  suggest  that 
the  theory  of  plasticity  be  used  to  represent  interfaces.  These  models  appear 
to  be  superior  to  the  conventional  Coulomb  model  for  several  reasons.  First, 
relative  motion  is  observed  for  any  tangential  force.  The  plasticity  model 
predicts  this  phenomenon;  the  Coulomb  model  does  not.  Second,  hardening  and 
softening  are  often  observed,  and  these  features  are  easily  characterized  by  a 
plasticity  model.  Third,  multidimensional  features  are  easily  incorporated  in 
a  plasticity  model,  whereas  special  rules  must  be  constructed  for  the  Coulomb 
model.  As  is  shown,  the  plasticity  formulation  can  also  be  incorporated  in 
existing  computer  codes  with  a  minimal  degree  of  adjustment. 

Ultimately  any  model  must  be  verified  by  experimental  data.  Some  data  on 
concrete-soil  interfaces  exist  (Ref.  25),  although  these  data  have  apparently 
never  been  used  for  verification  purposes.  Reference  25  gives  rather  exten¬ 
sive  results  for  both  clay  and  sand,  together  with  data  that  should  permit 
these  materials  to  be  well  defined. 
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The  cited  references  indicate  that  several  aspects  of  interface  modeling 
are  being  investigated  and  that,  as  usual,  a  large  disparity  exists  between 
important  theoretical  considerations  and  engineering  models. 
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V.  RECOMMENDATIONS 
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The  priniary  emphasis  In  this  work  phase  has  been  on  developing  an 
Improved  understanding  of  strain  invariants  and  rate  effects  and  on  Incorpora¬ 
ting  cracking  Into  the  constitutive  model.  Developments  associated  with 
strain  Invariants  have  been  quite  successful,  and  an  engineering  approach  for 
representing  cracking  has  been  demonstrated,  but  questions  about  rate  effects 
remain  unanswered.  The  basis  has  been  established  for  focusing  now  on  spe¬ 
cific  topics  that  address  fundamental  questions  concerning  material  behavior. 
The  answers  to  these  questions  could  provide  an  Improved  model.  The  topics  to 
be  addressed— rate  effects  in  concrete,  Interface  modeling,  and  the  response 
of  c1ay--have  been  selected  on  the  basis  of  the  general  dearth  of  knowledge  in 
these  areas. 

The  experimental  device  proposed  for  ascertaining  the  multiaxial  response 
of  concrete  under  strain  rates  having  a  magnitude  of  the  order  of  1/s  Is 
described  In  Section  II.  The  essential  aspect  of  the  device  is  that  the  com¬ 
bination  of  specimen  size,  strain  rate  magnitude,  and  loading  rates  obtained 
from  explosives  Is  uniquely  appropriate  In  that  all  major  criteria  can  be  met. 

The  experimental  equipment  can  help  answer  several  fundamental  questions. 
No  Information  is  available  on  whether  dilatation  Increases  or  decreases  with 
strain  rate.  The  enhancement  of  strength  with  strain  rate  as  a  function  of 
the  ratios  of  principal  stresses  Is  not  known.  With  a  plasticity  model ,  there 
Is  an  open  question  concerning  which  of  the  hardening  parameters—straln  rate. 
Inelastic  strain  rate,  or  stress  rate— Is  most  suitable  for  the  theoretical 
formulation.  Overall,  It  Is  believed  that  the  theoretical  basis  Is  suffi¬ 
ciently  advanced  to  provide  sound  guidelines  for  an  important  experimental 
program. 

The  observed  Interface  behavior  of  a  change  in  interface  shear  for  any 
change  In  relative  displacement  before  the  limit  state  is  reached  is  consider¬ 
ably  different  from  the  behavior  Implied  by  a  Coulomb  model.  The  use  of  the 
Coulomb  model  also  has  the  implication  of  possible  nonexistence  of  solutions. 
An  Interface  model  that  incorporates  newer  mathematical  concepts  such  as  a 
nonlocal  feature  may  be  especially  suitable  for  representing  the  roughness  of 
concrete  surfaces.  Furthermore,  such  an  approach  may  be  helpful  to  the  devel¬ 
opment  of  a  robust  numerical  algorithm. 
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The  viscoplastic  model  for  representing  frictional  materials  provides  a 
reasonable  basis  for  the  exploration  of  new  regimes  of  material  behavior.  I 
Is  recommended  that  the  model  be  used  as  a  basis  for  addressing  fundamental 
research  on  problems  concerning  rate  effects  in  concrete  and  soil -concrete 
interfaces. 
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APPENDIX  A* 


A  THIRD- INVAR I ANT  PLASTICITY  THEORY  FOR  LOW-STRENGTH  CONCRETE 


A  theory  of  plasticity  for  frictional  materials  is  developed  in  which 
first  and  third  invariants  of  stress  and  strain  are  used  instead  of  the  more 
conventional  second  invariants.  The  usual  concepts  of  strain-hardening  plas¬ 
ticity  are  used  with  the  exception  that  a  nonassociated  flow  rule  is  required 
to  control  dilatation.  A  procedure  for  determining  material  parameters  is 
outlined.  For  a  weak  concrete,  detailed  comparisons  are  made  between  theoret¬ 
ical  and  experimental  stress-strain  data  for  a  large  number  pf  three- 
dimensional  paths. 


♦This  appendix  is  a  reproduction  of  a  professional  paper  submitted  for  publi¬ 
cation  in  Journal  of  Engineering  Mechanics  (American  Society  of  Civil  Engi¬ 
neers).  Thus  it  is  a  self-contained  document  with  its  own  internally  consis¬ 
tent  numbering  system  for  equations,  references,  and  figures  and  with  a  format 
prescribed  by  the  publisher. 


A  THIRD- INVARIANT 


PLASTICITY  THEORY  FOR  LOW-STRENGTH  CONCRETE 

Howard  L.  Schreyer,1  M.  ASCE,  and  Susan  M.  Babcock,2  A.M.  ASCE 

September,  1983 


INTRODUCTION 

Numerous  attempts  have  been  made  to  model  concrete  with  the  theory  of 
plasticity.  Several  models  incorporate  strain  hardening,  but  none  has 
achieved  the  generality  and  simplicity  represented  by  the  Von  Mises  isotropic 
hardening  model  used  for  metals.  A  corresponding  theory,  in  which  the  first 
and  third  invariants  of  stress  and  strain  are  used  instead  of  the  second 
invariants,  is  developed  here  for  frictional  materials.  Material  parameters 
are  chosen  for  a  weak  concrete,  and  detailed  comparisons  are  made  between 
theoretical  and  experimental  stress-strain  data  for  a  large  nunber  of  three- 
dimensional  paths. 

In  the  phenomenological  approach  to  modeling  a  frictional  material  such 
as  concrete,  it  is  assuned  that  all  essential  response  features  can  be  repre¬ 
sented  in  a  tractable  model  without  the  explanation  that  would  be  provided  by 
micromechanics.  Such  an  approach  is  justified  by  the  need  to  provide  engi¬ 
neering  analyses  of  constructed  or  planned  systems.  The  resulting  model 
should  not  violate  basic  principles  and,  as  a  deeper  understanding  is 

Professor,  Department  of  Mechanical  Engineering,  and  Research  Scientist,  New 
Mexico  Engineering  Research  Institute,  University  of  New  Mexico, 

Albuquerque,  NM  87131. 

Research  Engineer,  New  Mexico  Engineering  Research  Institute,  University  of 
New  Mexico,  Albuquerque,  NM  87131. 


achieved,  additional  details  should  be  incorporated  easily  within  the  frame¬ 
work.  However,  the  current  status  of  model  development  requires  that  the  ade 
quacy  of  a  model  be  evaluated,  not  on  a  micromechanical  basis,  but  on  a  macro 
mechanical  basis  as  represented  by  mechanical  experiments  on  specimens  in  a 
laboratory.  For  concrete  loaded  in  compression  to  the  limit  state,  the  gen¬ 
eral  features  exhibited  include  inelastic  compaction  under  hydrostatic  load¬ 
ing,  shear  enhanced  compaction,  dilatation,  and  strain  hardening.  These 
aspects  of  stress-strain  behavior  must  be  simulated  realistically  to  obtain 
engineering  solutions  to  a  variety  of  boundary-value  problems  involving  con¬ 
crete  structures. 

i 

Representative  models  for  concrete  are  described  in  Refs.  1,  2,  3,  and  8 
Often  the  primary  concern  is  the  limit  state,  in  which  case  an  elastic-per- 
fectly  plastic  approach  is  adequate  (8).  In  other  cases,  the  actual  strain¬ 
hardening  features  are  represented.  Comparison  with  experimental  data  are 
often  provided,  but  for  the  most  part  these  data  are  of  a  restricted  class, 
triaxial  paths  or  biaxial  paths  in  plane  stress  being  the  most  common.  As 
such,  these  data  can  provide  only  a  limited  evaluation  of  models. 

In  a  recent  set  of  experiments  on  cubical  specimens,  three-dimensional 
paths  were  prescribed  (6).  The  paths  varied  from  triaxial  to  multi  segmental 
paths  that  provided  transitions  between  uniaxial,  biaxial,  and  proportional 
loadings.  Each  path  was  performed  on  two  or  three  specimens.  The  results 
were  evaluated  and  combined  into  single  stress-strain  curves  for  each  path 
(7).  Although  principal  axes  of  stress  were  held  fixed,  it  is  believed  that 
these  data  can  provide  a  more  comprehensive  evaluation  of  a  model  than  that 
normally  undertaken. 


In  Ref.  4  a  third-invariant  model  for  frictional  materals  is  described, 
and  an  initial  evaluation  based  on  a  comparison  of  theoretical  and  experimen¬ 
tal  data  for  two  materials  is  provided.  Analogous  to  the  von  Mises  hardening 
theory,  the  model  is  an  isotropic  strain-hardening  one;  therefore,  there  are 
relatively  few  parameters.  Thus,  if  the  basic  response  features  of  concrete 
can  be  adequately  represented,  the  model  would  be  appropriate  for  use  in  large 
finite-element  or  finite-difference  codes.  An  evaluation  is  provided  by  using 
a  large  number  of  paths  for  one  material,  rather  than  a  few  paths  for  several 
materials.  In  addition  to  the  correlation  between  theoretical  and  experimen¬ 
tal  results,  the  procedure  for  choosing  material  parameters  is  described. 

EXPERIMENTAL  DATA 

The  triaxial  testing  device  at  New  Mexico  State  University  can  apply, 
independently,  loads  on  three  faces  of  a  cube  so  that  all  paths  are  stress- 
controlled.  This  device  was  used  to  obtain  stress-strain  data  for  a  variety 
of  paths  (6).  For  each  path  these  data  were  evaluated  and  condensed  to  repre¬ 
sentative  curves  (7),  which  can  be  used  conveniently  for  evaluating  constitu¬ 
tive  models. 

The  tests  were  divided  into  three  sets,  the  first  two  of  which  were  con¬ 
structed  with  all  principal  components  of  stress  in  compression.  The  last 
set,  which  will  not  be  considered  in  this  paper,  involved  a  tensile  component 
of  stress  for  at  least  some  portion  of  each  path.  Each  test  in  the  first  set 
was  performed  on  two  specimens.  The  tests  consisted  primarily  of  conventional 
uniaxial,  biaxial,  triaxial,  and  shear  tests.  The  second  set  of  nine  tests 
involved  more  complicated  paths  In  the  compressive  regime;  consequently,  each 
test  in  this  set  was  performed  on  three  specimens. 


To  perform  each  experiment,  it  was  necessary  to  give  the  stress  paths  in 
incremental  form,  which  is  also  a  convenient  way  to  describe  the  paths  analy¬ 
tically.  Tables  1  and  2  provide  a  complete  description  of  these  paths  in 
terms  of  a  positive  increment  in  stress,  As.  Each  path  is  identified  by  a 
number  and  the  set  to  which  it  belongs,  and  each  segment  of  the  path  is 
described  briefly.  Terms  are  defined  in  the  tables  as  necessary. 

The  paths  are  shown  in  Fig.  1  for  the  P  plane  where 

P  =  -  ^(o!  +  o2  +  o3 )  ) 

^  ■{£(<>,  -°2>2  *  (oj  -»3)2  +  («3  "  Oi)2]}1/2  (2) 

represent  the  mean  pressure  and  the  second  invariant  of  the  stress  deviator. 
The  principal  stresses,  op  o2»  and  o3,  are  assuned  positive  in  tension.  As 
these  plots  vividly  demonstrate,  the  experiments  cover  a  large  spectrum  of 
stress -defined  paths. 

Considering  the  complexity  of  the  paths,  a  remarkable  repeatability  was 
obtained.  Certain  anomolous  results  were  observed,  however,  and  corrections 
were  made  for  obvious  shifts  in  strain  data.  All  such  adjustments  are  fully 
described  in  Ref.  7.  After  these  adjustments  were  made,  composite  stress- 
strain  curves  were  developed  for  each  stress  path.  These  data  are  given  as 
the  experimental  curves  shown  in  Fig.  2.  To  illustrate  the  large  variations 
in  stress  and  strain  from  one  path  to  the  next,  the  same  scale  is  used  for  all 
the  plots.  To  conserve  space,  any  curve  starting  from  the  lower  right  corner 
is  associated  with  the  abscissa  at  the  top  of  each  plot.  The  variation  of 
pressure  with  volumetric  strain  is  a  useful  relation,  but  this  information  is 
implicit  in  Fig.  2. 


TABLE  1.  SET  1:  STRESS  PATHS  IN  COMPRESSION 


Path 

1 


Uniaxial  stress  to  failure 
A®!  3  Ao2  *  0;  Ao3  *  -As 


Description 


2 


Uniaxial  stress  to  failure.  This  and  path  1  represent  a  test 
for  initial  anisotropy. 
aol  *  -as;  ao2  3  a<j3  3  0 


3 


5 


6 


7 


8 


Biaxial  stress  to  failure  with  o2/oi  =  1.0;  c3  3  0 

Aoi  3  -as;  a o2  *  -as;  ao3  *  0 
% 

Hydrostatic  loading  followed  by  uniaxial  loading  to  failure 
(1)  Aci  3  Ao2  3  A03  3  -As 

to  the  point  0l  3  o2  *  o3  =  -8.5  MPa 
(i 1 )  Aa:  3  -as;  a o2  3  ao3  *  0 

Hydrostatic  loading  followed  by  shear  loading  to  failure 
(i)  Aoa  3  ha 2  3  Ao3  3  “AS 

to  the  point  ox  3  o2  3  a3  3  -16.6  MPa 
(ii)  Aoi  3  "as;  ao2  -  as;  ao3  -  0 
3-D  stress  to  failure  (subject  to  extensometer  capability)  with 
a2/o!3  1.0;  o3 loi  3  0.1 
Aai  3  -as;  ao2  3  -as;  Aa3  3  -0.1  As 
Biaxial  stress  to  failure  with  o2/oi  3  0.5;  o3  3  0 
Ao}  3  -as;  ao2  3  -0.5  as;  a<j3  3  0 

3-0  stress  to  failure  (subject  to  extensometer  capability)  with 

o2/oi  3  0.5;  03/Ol  3  0.1 

Aoi  3  -as;  Ao2  3  -0.5  as;  Aa3  3  -0.1  As 


TABLE  2.  SET  2:  ADDITIONAL  STRESS  PATHS  IN  COMPRESSION 


Description 


(i) 

Load 

hydrostatically  to  o:  -  a2  - 

Q 

CO 

11 

1 

2.2  MPa 

bo±  =  bo  2 

3  ao3 

3  -AS 

(ID 

Load 

uni  axially  to  oj.  *  -17.8  MPa 

;  o2  = 

o3  3  -2.2 

MPa 

Aa* 

3  -as;  a©2  =  A03  -  0 

(HD 

Unload  to  <?!  s  a2  s  a 3  =  “2.2  MPa; 

Ao  1  3 

as;  a 02  3 

Ao3  3  0 

(iv) 

Load 

biaxially  in  xi-x2  plane  to  failure 

with 

a2/a 

l  3  1.0;  a3  3  -2.2  MPa 

Aa* 

=  A02  ~  -as;  A03  -  0 

(1) 

Load 

hydrostatically  to  0l  =  a2  «  , 

c3  3  - 

2.1  MPa; 

=  A02  “  Aa 3  *  -As 

(ID 

Load 

biaxially  in  xz-x2  plane  to  failure 

with  al/a2  ~  1*0; 

°3  * 

-2.1  MPa 

Ao^ 

*  Ao2  *  -as;  act3  3  0 

(1) 

Load 

hydrostatically  to  Oi  3  a2  3 

a3  3  - 

2.0  MPa 

AaA 

3  A a2  *Ao3  3  -AS 

(H) 

Biaxial  stress  in  x1-x2  plane  to 

0 1  * 

o2  *  -22.5  MPa;  0 3  3  -2.0  MPa 

Aoi 

3  Ao2  *  -as;  ao3  *  0 

(111) 

Unload  to  aj  *  o2  *  03  3  -2.0  MPa; 

Aai  3 

A a2  3  as; 

Ao3  *  0 

(Iv) 

Load 

uniaxially  in  x: -direction  to 

failure 

Aoi 

3  -as;  Aa2  3  A03  3  0 

(D 

Load 

hydrostatically  to  Oi  =02  = 

03  3  " 

2.0  MPa 

Aai 

3  Ao2  3  Aa3  *  -As 

(ID 

Load 

uniaxially  In  Xj -direction  to 

failure 

A  0! 

3  as;  ao2  =  Ao3  3  0 

TABLE  2.  CONTINUED 


Path  Description 

5  (i)  Load  hydrostatically  to  <j1  s  a2  a  03  3  -26.2  MPa 

Acj  *  Ao2  *  A03  *  -As 

(ii)  Unload  hydrostatically  to  ^  3  o2  3  o3  3  "17.5 

A01  *  A<J2  *  A03  a  4s 
(iii)  Load  in  shear  to  failure 

AO}  3  —AS}  ao2  *  AS}  Adj  =  0 

6  (i)  Load  hydrostatically  to  -  a2  ~  03  3  -40.0  MPa 

Aoa  3  Ao2  3  A03  a  -As 

(ii)  Unload  hydrostatically  to  s  a2  8  a3  s  -2.0  MPa 
A<jj  3  A 02  3  AO3  3  As 
(iii)  Load  biaxial ly  to  failure 


Aot  a  Ao2  3  -AS}  Ao3  3  0 


(1) 

Load 

biaxial ly  to  ot  3  o2  3 

-2.2  1 

MPa;  a3 

3  0 

Aol  = 

s  Ao2  3  -AS}  ao3  3  0 

(Ii) 

Load 

uniaxially  to  Oi  3  -15.1 

3  MPa 

;  a2  a 

-2.2  MPa} 

o3  3  0 

AOi  3 

*  -AS}  ao2  a  A03  3  0 

(iii) 

Load 

to  a  biaxial  state  a 

°2  S 

-15.8 

MPa}  o3  * 

0 

AOj 

5  a<73  *  0}  a o2  *  -4s 

(iv) 

Load 

biaxially  to  failure 

Aoj  3 

s  Ao2  3  -as;  Aa3  3  0 

(i) 

Load 

hydrostatically  to  3 

02  3 

®3  =  - 

2.1  MPa 

Aoj  3  A a2  a  A03  3  -As 


(ii)  Load  uniaxially  to  ol  3  -17.6  MPa 
Aox  3  -AS}  Ao2  3  a<j3  *  0 

(iii)  Trace  a  path  to  3  -2.1  MPa}  o2  3  -17.6  MPa}  o3  3  -2.1  MPa 
A01  3  AS}  Ao2  *  -AS}  Ao3  3  0 


TABLE  2.  CONCLUDED 


Description 


(iv) 

Load 

uni  axially  to 

failure 

Aoi 

~  ha 3  =  0;  A02 

=  -As 

(i) 

Load 

hydrostatically  to  0l 

Aci 

“  ha 2  ~  Ao 3  “ 

-As 

(ii) 

One 

cycle  in  shear 

Stress,  MPa  Stress,  MPa 


Stress,  MPt 


Stress,  MPa 


Figure  2,  Concluded. 


THIRD- INVAR I ANT  MODEL 


The  third-invariant  model  has  been  described  in  Ref.  4,  but  the  theory 
has  been  modified  slightly  since  then.  Thus,  the  formulation  is  outlined  here 
but  with  more  emphasis  placed  on  the  evaluation  of  material  parameters. 

A  particular  combination  of  the  three  stress  invariants  can  also  be  con¬ 
sidered  a  third  invariant.  On  the  basis  of  the  use  of  a  shifted  stress  con¬ 
cept,  such  an  invariant  is  defined  by 

L  =  -  [(ax  -  cs)(a2  -  crs ) (ct3  -  o$)  +  a|]1/3  (3) 

in  which  the  parameter  o$  can  be  defined  as  the  cohesion  and  is  assumed  to 
vary  with  the  deformation.  On  the  basis  of  previous  fits  to  experimental 
data,  the  limit  surface  is  assumed  to  satisfy  the  relations 

°s  *  °L  (4> 

L  ■  YLp  -  <*L  (5) 

Because  limit  point  data  for  a  set  of  paths  are  available,  yl  and  cr^  can  be 
obtained  such  that  the  limit  line  is  an  optimal  fit  in  a  least  squares  sense. 
For  consistency,  the  value  of  must  be  chosen  so  that  the  value  used  in 
Eq.  3  agrees  with  the  intercept  value  of  the  straight  line  defined  by  Eq.  5. 
Then  yl  is  chosen  to  maximize  the  coefficient  of  correlation.  The  result  of 
using  the  data  of  sets  1  and  2  is  aL  -  1.87  MPa  and  yl  =  0.79,  with  a  coeffi¬ 
cient  of  correlation  of  0.981.  The  corresponding  fit  is  shown  in  Fig.  3. 

The  flow  function  is  defined  by 

♦  ■  YP  -  L  -  a$  (6) 

so  that  the  flow  surface  $  *  0  is  also  a  line  in  the  L-P  space.  The  strain¬ 
hardening  functions  y  and  a  are  assumed  to  be  similar  in  form.  As  hardening 


Pressure  (P),  MPa 


Figure  3.  Straight-line  fit  to  limit-point  data. 

occurs,  y  decreases  from  an  initial  value  of  yQ  to  the  limit  value  y^,  while 
0  increases  from  o0  to  o^.  Softening  is  achieved  by  allowing  y  to  increase 
and  o$  to  decrease  after  the  limit  state.  The  predicted  response  is  not 

overly  sensitive  to  the  initial  values,  which  are  obtained  by  noting  when 
nonlinearity  first  appears  on  a  typical  stress-strain  path.  These  values  were 

chosen  to  be  yo  a  2.5  and  o0  *  1.0  MPa. 

To  determine  suitable  functions  for  y  and  os,  suppose  that  a  path  length 

strain  invariant,  e**,  exists  and  that  it  has  a  unique  value,  e^,  at  the  limit 
state  no  matter  which  path  is  followed.  Appropriate  strain-hardening 
functions  that  automatically  provide  the  correct  initial  values  of  t  *  y0  and 


o$  s  ag  when  e  =  0  and  limit  state  values  of  y  =  yl  and  a$  =  when 


T1  =  e|[  are  then  given  by 


Y  *  Y0  +  (yL  -  Y0)  sin  U 

2W. 


(7) 


it  /  e1  ' n 

os  *  O0  +  (°L  -  o0)  sin  -  — 

2W 


(8) 


The  parameter  n  is  chosen  to  provide  a  reasonable  shape  for  the  strain¬ 
hardening  portion  of  the  predicted  stress-strain  curves.  A  value  of  n  *  0.16 
was  considered  appropriate  for  the  weak  concrete  on  the  basis  of  comparisons 
between  theoretical  and  experimental  stress-strain  curves. 

Because  an  associated  flow  rule  provides  too  much  dilatation,  a  potential 
function  of  the  following  form  is  used: 

<J>*  *  y*P  -  L  -  o  (9) 

where 


Y*  3  y  +  b  (10) 

and  b  is  a  material  constant.  A  nonassociated  flow  rule  implies  that  the 
inelastic  strain  increment  is  given  by 


01) 


in  which  dx  is  a  positive  parameter,  and  a  and  I  denote  the  stress  and  iden¬ 
tity  tensors,  respectively.  The  implication  of  this  equation  is  that  dilatation 


is  reduced  from  that  given  by  the  associated  flow  rule  if  b  is  positive.  For 
hydrostatic  compaction  with  P  »  a ,  it  follows  that  L  *  P,  and  the  volumetric 
inelastic  strain  increment  is  given  by 

de’  *  trace  dg1 

*  dX(l  -  y-  b)  (12) 

The  requirement  that  such  a  loading  path  should  always  yield  compaction  is 
satisfied  if  b  >  1  -  Y|_.  The  easiest  way  to  determine  b  is  to  monitor  the 
lateral  strains  on  a  typical  triaxial  loading  path  and  to  adjust  b  to  provide 
the  type  of  fit  desired  for  an  ensemble  of  such  paths.  Here,  the  value  of 
b  =  0.35  was  obtained. 

An  appropriate  form  for  F1  is  still  open  to  question,  but  such  a  form 
must  include  a  volumetric  effect  if  inelasticity  is  to  be  obtained  for  hydro¬ 
static  loading.  Both  compaction  and  dilatation  can  occur;  therefore,  an 
absolute  value  is  used  to  obtain  the  path  length  increment  of  the  first 
invariant: 

d¥{  *  tr  de^  (13) 

Because  the  second-invariant  increment  is  so  useful  for  the  theory  of 
metal  plasticity,  it  is  included  as  follows: 

del  ■  [tr(d|,d)2]1/2  (14) 

In  which  de™  denotes  an  Increment  in  the  inelastic  strain  deviator. 


For  the  sake  of  generality  a  third  invariant  increment  is  also  intro¬ 


duced,  with  the  absolute  value  used  to  ensure  a  monotonically  increasing 
parameter: 

die]  *  |  detfde1)  |  1/3  (15) 

An  increment  in  a  path  length  invariant  for  frictional  materials  is  assumed  to 
be 


ci  del  +  C2del  +  C3de] 


(16) 


where  clt  c2,  and  c3  are  material  parameters.  For  any  path,  it  is  possible  to 
determine  the  value  of  e  at  the  limit  point;  the  limit  value  is  denoted  by 
e’*.  If  the  invariant  e1*  is  to  provide  a  useful  plasticity  theory,  then  "e** 
should  be  the  same  for  all  paths.  Data  such  as  those  provided  in  Refs.  6  and 

7  can  be  used  to  determine  whether  or  not  "ej*  is  unique  for  a  given  material. 

In  fact,  when  these  data  are  used,  values  for  e^*vary  by  an  order  of  magnitude 
no  matter  which  values  are  chosen  for  cx,  c2,  and  c3.  However,  a  distinguish- 

— *i  * 

ing  feature  does  appear:  for  those  paths  associated  with  large  P,  e  is  much 
larger  than  the  value  obtained  for  uniaxial  stress.  This  dependence  on  P 
suggests  that  the  definition  of  the  third  invariant  should  involve  a  weighting 
function.  On  the  basis  of  a  preliminary  evaluation  (5),  the  following  form 
for  an  increment  is  proposed: 


Po 

P  +  P0 


(17) 


Here,  the  constant  P0  controls  the  effect  of  mean  stress.  If  P0  is  large,  the 
effect  of  P  is  minimal;  the  reverse  is  true  if  P  is  small. 


The  use  of  e1  rather  than  e1*  produces  a  dramatic  reduction  in  the  varia¬ 
tion  of  values  for  "e*  at  the  limit  state,  suggesting  that  the  parameters  c*, 

Co ,  c,,  and  Pn  can  be  adjusted  to  minimize  this  variation.  Because  the  nature 


of  the  data  hardly  warrants  a  formal  minimization  procedure,  a  simple  search 
procedure  was  instigated  in  which  one  parameter  was  adjusted  at  a  time.  With 
the  arbitrary  choice  of  ci  =  1,  the  resulting  values  for  the  other  parameters 
are  P0  =  7  MPa,  c2  =  0,  and  c3  =  0.5.  The  interesting  implication  of  c2  =  0 
is  that  only  the  first  and  third  invariants  of  strain  are  necessary,  which 
provides  a  nice  correlation  with  the  use  of  the  first  and  third  invariants  of 
stress. 

The  resulting  values  for  at  the  limit  states  for  all  paths  are  given 
in  Table  3.  At  first  glance  the  variation  still  appears  to  be  exceedingly 
large,  although  if  one  or  two  anomalous  points  such  as  paths  3  and  7  of  set  2 

are  excluded,  this  variation  is  reduced  considerably.  Paths  1  and  2  of  set  1 

are  both  uniaxial  stress  paths,  which  were  performed  to  check  anisotropy  by 
testing  along  different  axes  on  two  specimens  from  the  same  batch.  The  corre¬ 
sponding  values  of  "e^  differ  by  a  factor  of  2.  It  was  also  noted  by  Traina 
(6)  that  at  the  limit  points,  the  strains  were  increasing  at  a  significant 

rate  and  that  it  was  difficult  to  identify  closely  the  strains  at  the  limit 

point.  Furthermore,  the  data  were  obtained  directly  in  digital  form  at  each 
of  approximately  ten  steps  for  each  path.  This  digitizing  can  introduce  sig¬ 
nificant  error  through  the  numerical  sum  used  to  approximate  the  integral  for 

e1.  For  these  reasons,  more  sophisticated  approaches  to  evaluating  Po  and  C3 
from  these  data  are  unwarranted.  Furthermore,  it  appears  that  predicted 

results  are  not  overly  sensitive  to  e^;  therefore,  a  representative  value  of 

e£  s  0.0025  was  selected. 

On  the  basis  of  a  combination  of  classical  approaches  and  initial  loading 
characteristics,  the  remaining  parameters  of  Young's  modulus  and  Poisson's 
ratio  were  chosen  to  be  E  *  21,000  MPa  and  v  *  0.25,  respectively. 


TABLE  3.  VALUES  OF  THE  STRAIN  INVARIANT  V  AT  LIMIT  POINTS 


RESULTS 

With  the  material  parameters  identified  in  the  previous  section,  theoret¬ 
ical  predictions  of  strain  were  made  for  the  stress  paths  listed  in  Table  1. 
These  results  are  shown  with  the  experimental  data  in  Fig.  2.  In  general,  the 
correlation  is  good  considering  that  the  model  is  that  of  isotropic  strain¬ 
hardening  plasticity. 

Variations  between  theoretical  and  experimental  values  of  peak  stress  are 
a  direct  consequence  of  how  well  the  straight  line  of  Fig.  3  fits  the  limit 
state  data.  Corresponding  variations  in  strain  at  the  limit  state  are  reflec¬ 
ted  to  a  considerable  extent  in  how  well  the  invariant  i"1  captures  the  physi¬ 
cal  phenomena.  Strains  at  the  limit  state  are  also  affected  by  the  path; 
therefore,  agreement  is  contingent  on  a  faithful  representation  through  the 
strain-hardening  portion  of  the  path  history. 


The  largest  discrepancy  occurs  for  the  proportional  loading  case  given  as 
path  8  of  set  1.  The  predicted  peak  stress  is  much  higher  than  the  experimental 
value  even  though  the  limit  condition  is  satisfied.  The  reason  is  that  in  the 
L-P  plane,  a  proportional  path  such  as  this  one  is  almost  parallel  to  the 
limit  line  so  that  a  slight  variation  is  magnified  considerably.  If  the  limit 
line  were  curved  rather  than  straight,  the  discrepancy  might  be  removed. 

Also,  the  measured  intermediate  strain  is  positive  (tensile),  whereas  the 
predicted  strain  is  negative.  A  similar  feature  is  displayed  to  a  lesser 
extent  for  path  5,  sets  1  and  2.  The  reason  is  that  the  shape  of  the  flow 
surface  in  the  pi-plane,  as  defined  implicitly  through  the  flow  condition 
♦  *  0,  is  not  as  accurate  as  that  proposed  by  Wiliam  and  Warnke  (8).  In  other 
words,  an  error  is  introduced  when  a  simplified  model  is  used.  Whether  or  not 
the  error  is  significant  depends  on  the  intended  application.  To  illustrate 
the  point  further,  suppose  path  8  of  set  1  is  altered  so  that  the  intermediate 
stress  is  0.25  times  the  primary  stress.  This  alteration  corresponds  to  a 
slightly  different  radial  line  in  the  pi-plane.  The  model  predicts  essen¬ 
tially  zero  intermediate  strain  as  shown  in  Fig.  4  so  that  the  theoretical 
ratio  of  0.25  is  the  transition  value  between  positive  and  negative  interme¬ 
diate  strains.  The  higher  value  indicated  by  the  experimental  data  may  be  the 
result  of  small  shear  stresses  on  the  loaded  faces  of  the  specimens. 

The  next  most  important  discrepancy  involves  the  amount  of  inelastic 
deformation  under  hydrostatic  loading  as  shown,  for  example,  in  the  first  part 
of  path  9,  set  2.  The  predicted  inelastic  volumetric  strain  is  less  than  the 
experimental  strain.  Also,  upon  unloading,  the  model  always  predicts  an  elas¬ 
tic  response,  whereas  some  inelastic  deformation  is  displayed  in  the  tests.  A 
kinematic  hardening  aspect  would  have  to  be  added  to  the  model  to  enable  it  to 
predict  such  a  response  feature. 
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Figure  4.  Theoretical  prediction  for  triaxial  loading 
with  o2  =  0.25crt  and  *  O.loj. 


For  all  other  cases  it  is  believed  that  the  agreement  is  more  than  ade¬ 


quate  for  the  purposes  of  most  engineering  analyses. 


CONCLUSION 

These  results  show  that  there  is  considerable  merit  in  the  use  of  a 
strain-hardening  isotropic  plasticity  model  for  representing  the  essential 
features  of  concrete  behavior  in  the  compressive  regime.  The  study  indicates 
that  the  first  and  third  invariants  of  stress  and  strain  more  appropriately 
represent  concrete  than  do  the  second  invariants,  which  are  commonly  used  for 
metal  plasticity.  A  nonassociated  flow  rule  is  necessary  to  control 
dilatation. 


Several  aspects  of  concrete  behavior  have  not  been  considered.  Among 
these  are  tensile  behavior,  strain  softening,  and  rate  effects.  The  current 


theory  has  been  extended  to  incorporate  these  effects,  but  until  more  multi- 
axial  experimental  data  become  available,  no  definitive  evaluation  can  be 
made.  The  primary  advantage  of  the  extended  model  is  that  it  provides  good 
guidelines  for  developing  experimental  apparatus  with  appropriate 
instrumentation. 

Because  the  formulation  is  patterned  after  correspondingly  attractive 
models  used  for  metals,  a  constitutive  equation  algorithm  can  easily  be  added 
to  existing  computer  codes.  The  model  implicitly  incorporates  features  such 
as  inelastic  hydrostatic  compaction,  shear  enhanced  compaction,  dilatation, 
and  strain-hardening.  For  those  engineering  applications  in  which  a  predic¬ 
tion  of  deformation  is  essential,  this  theory  provides  a  balanced  approach 
between  accuracy  and  computational  feasibility.  Although  the  results  of  this 
work  are  restricted  to  low  strength  concrete,  the  model  should  be  practicable 
for  other  concretes  and  geological  materials. 
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NOTATION 

b  -  material  parameter; 


Cj,  c2.  c3  *  material  parameters; 
dA  =  material  parameter; 


E  =  Young's  modulus; 
e3  =  principal  strains; 

e^  =  inelastic  volumetric  strain; 

e°  =  path  length  strain  invariant; 

"e1*  =  path  length  strain  invariant; 
e^  =  path  length  strain  invariant  limit  value; 
g1  =  inelastic  strain  tensor; 

gld  =  inelastic  deviatoric  strain  tensor; 

J  =  identity  tensor; 

I 

J2  =  second  invariant  of  the  stress  deviator; 

L  *  stress  invariant; 
n  =  material  parameter; 

P  =  mean  pressure; 

P0  =  material  parameter  with  units  of  pressure; 
Y  -  slope  of  flow  surface  in  L-P  space; 
yl  =  slope  of  limit  surface  in  L-P  space; 

Y0  =  initial  slope  of  flow  surface  in  L-P  space 
Y*  =  slope  of  potential  function  in  L-P  space; 
As  =  positive  increment  in  stress; 
v  =  Poisson's  ratio; 

=  cohesion  at  limit  surface; 

*  cohesion; 

o0  3  initial  cohesion; 
o 3  =  principal  stresses; 

0  =  stress  tensor; 

$  =  flow  function;  and 
a  potential  function. 


APPENDIX  B* 


A  THIRD- INVAR  I ANT  VISCOPLASTICITY  THEORY  FOR  RATE-DEPENDENT  SOILS 


A  theory  of  viscoplasticity  for  frictional  materials  is  developed  in 
which  first  and  third  invariants  of  stress  and  strain  are  used  instead  of  the 
more  conventional  second  invariants.  Rate  effects  are  incorporated  directly 
into  the  expression  for  the  flow  surface  so  that  the  numerical  algorithm  for  a 
plasticity  subroutine  can  be  used.  The  usual  concepts  of  strain-  and  strain 
rate-hardening  viscoplasticity  are  used  except  that  a  nonassociated  flow  rule 
is  required  to  control  dilatation.  For  two  sandy  materials,  detailed  compar¬ 
isons  are  made  of  theoretical  and  experimental  stress-strain  data  for  both 
static  and  dynamic  paths. 


♦This  appendix  is  a  reproduction  of  a  professional  paper  submitted  for  publi¬ 
cation  in  Journal  of  Engineering  Mechanics  (American  Society  of  Civil  Engi¬ 
neers).  Thus  it  is  a  self-contained  document  with  its  own  internally  consis¬ 
tent  numbering  system  for  equations,  references,  and  figures  and  with  a  format 
prescribed  by  the  publisher. 


A  THIRD-INVARIANT  VISCOPLASTICITY  THEORY  FOR 
RATE-DEPENDENT  SOILS 

Howard  L.  Schreyer,1  M.  ASCE,  and  James  E.  Bean,2  M.  ASCE 

October  1983 

INTRODUCTION 

Large  strain  rates  in  soils  arise  in  a  number  of  situations  that  include 
strong  seismic  motion,  pile  driving,  and  high-explosive  testing.  Displace¬ 
ments  as  a  function  of  time  for  points  in  the  earth  close  to  these  disturb¬ 
ances  are  frequently  required,  and  in  these  cases  the  constitutive  model  of 
the  ground  media  must  include  rate  effects.  With  modern  computer  codes, 
three-dimensional  analyses  are  feasible.  Thus,  the  model  must  be  quite  gen¬ 
eral  but  at  the  same  time  sufficiently  simple  so  that  storage  and  computation 
time  requirements  are  reasonable  for  numerical  computations.  Proposed  in  this 
paper  is  a  three-dimensional  engineering  model  that  provides  a  balance  between 
accuracy  and  computational  efficiency.  Material  parameters  are  selected,  and 
detailed  comparisons  between  theoretical  and  experimental  stress-strain  data 
are  given  for  a  sand  and  a  clayey  sand  subjected  to  both  static  and  dynamic 
loads. 

The  rate  effect  known  as  creep  involves  a  relatively  large  time  scale. 

The  rate  effect  of  interest  here  occurs  over  a  short  time  and  is  characterized 
by  an  apparent  enhancement  of  soil  stiffness  and  strength.  Materials  exhibit 
both  types  of  rate  effects,  although  the  characteristic  times  can  vary  widely. 
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Consequently,  statements  have  been  made  to  the  effect  that  particular  materi¬ 
als,  such  as  sandy  soils,  are  rate-independent.  However,  when  the  strain  rate 
is  large  enough,  a  significant  change  from  static  behavior  is  observed  even  in 
these  materials. 

Jackson  et  al.  (4)  at  the  Waterways  Experiment  Station  (WES)  performed  a 
set  of  static  and  dynamic  tests  on  a  clayey  sand.  The  three  principal  compo¬ 
nents  of  stress  and  strain  were  measured  in  the  static  tests,  and  the  result¬ 
ing  data  are  invaluable  for  the  assessment  of  constitutive  models.  In  the 
dynamic  tests,  which  were  uniaxial  strain  tests,  the  lateral  stresses  were  not 
measured;  therefore,  these  data  are  not  as  complete  as  desired.  The  fact 
remains,  however,  that  these  data  are  unique  because  so  little  is  available  in 
the  way  of  experimental  data.  An  exception  is  a  paper  by  Ito  and  Fujimoto  (3) 
in  which  experimental  results  are  given  for  sand  subjected  to  triaxial  path 
tests.  Because  characteristic  strain  rates  for  these  two  test  programs  vary 
from  1/s  (3)  to  200/s  (4),  the  data  from  these  experiments  provide  a  severe 
test  for  a  constitutive  model. 

Various  cap  models  have  evolved  from  the  original  work  of  Drucker  et  al. 
(2),  and  considering  the  simplicity  of  the  cap  formulation,  predictions  are 
remarkably  good.  However,  as  the  three-dimensional  response  features  of  soils 
have  become  better  understood,  the  models  have  had  to  be  modified.  Nelson  et 
al.  (7)  showed  how  a  hardening  law  could  be  established  to  control  dilatation, 
among  other  features.  This  approach  has  been  extended  by  Baladi  and  Rohani 
(1)  to  a  viscoplastic  form  that  represents  a  generalization  of  a  theory  given 
by  Perzyna  (8) . 
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The  cap  models  are  generally  formulated  in  terms  of  the  first  and  second 
invariants  of  stress  and  strain*  which  implies  that  the  shape  of  the  limit 
surface  in  the  pi-plane  is  circular.  However,  it  is  believed  (6)  that  for 
geological  materials  the  limit  surface  is  triangular  rather  than  circular; 
therefore,  a  more  accurate  theory  necessary,  and  the  third  invariant  of  stress 
must  be  used  in  addition  to  the  other  two.  This  complication  has  been  obvi¬ 
ated  to  a  certain  extent  by  the  observation  of  Lade  and  Duncan  (6)  that 
equally  adequate  predictions  of  limit  states  can  be  obtained  by  using  only  the 
first  and  third  invariants.  This  idea  has  been  developed  by  Lade  (5)  and  more 
recently  by  Schreyer  (11),  who  have  used  various  forms  of  the  third  invariant. 

Schreyer  and  Babcock  (12)  have  extended  the  concept  to  a  more  conventional 

strain-hardening  formulation  in  which  a  single  path-length  strain  invariant  is  1 

shown  to  represent  adequately  a  variety  of  paths  for  a  weak  concrete.  Their 

basic  formulation  is  shown  here  to  be  suitable  for  sandy  soils  as  demonstrated 

by  comparisons  of  theoretical  and  experimental  data.  In  addition,  the  theory 

has  been  extended  to  viscoplasticity,  but  by  means  of  an  approach  that  differs 

considerably  from  that  used  by  Baladi  and  Rohani  (1).  Instead  of  using  a 

viscosity  parameter,  the  approach  outlined  here  involves  describing  the  flow 

surface  in  terms  of  total  strain  rate  (9).  The  use  of  the  total  strain  rate 

provides  an  enhanced  region  of  elastic  response.  The  sets  of  experimental 

data  that  involve  rate  effects  are  used  to  evaluate  the  model. 

THIRD-INVARIANT  MODEL 

The  model  described  in  Ref.  11  has  been  modified  to  incorporate  rate 
effects.  For  the  sake  of  completeness,  the  basic  aspects  of  the  model,  with 
special  emphasis  on  sandy  soils,  are  summarized  in  this  section. 

I 
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Any  combination  of  three  independent  stress  invariants  is  also  an  invar¬ 
iant.  On  the  basis  of  a  shifted  stress,  first  and  third  stress  invariants  are 
defined  as  follows: 

P  =  “T  (°i  +  a2  +  a3 )  (!) 

L  *  -  [(oi  -  as)(o2  -  o$)(a3  -  ag)  +  a|]^3  (2) 

in  which  ox,  a2,  and  a3  are  principal  components  of  stress,  and  is  a  param¬ 
eter  that  can  be  defined  as  the  cohesion.  For  physical  reasons  and  also  for 
% 

the  sake  of  numerical  stability,  is  taken  to  be  a  small  positive  stress. 

On  the  basis  of  previous  results  (11),  the  limit  surface  (line)  is  assumed  to 
satisfy  the  relation 

L  *  7LP  -  o$  (3) 

in  which  yl  is  obtained  so  that  the  limit  line  is  an  optimal  fit  to  the  data 
in  a  least-squares  sense. 

A  flow  function  is  chosen  to  be 

<t>  =  YP  -  L  -  o$  -  (4) 

in  which  y  is  a  strain-hardening  function.  The  additional  function,  o  , 
depends  only  on  strain  rate  and  can  be  interpreted  as  an  enhancement  of  cohe¬ 
sion  based  on  rate  effects.  The  theory  falls  within  the  general  category 
considered  by  Rubin  (9). 
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The  flow  surface  ♦  5  0  is  also  a  line  in  the  L-P  space.  As  hardening 
occurs,  y  decreases  from  an  initial  value  of  Yo  to  a  limit  value,  y^.  Soften¬ 
ing  is  modeled  by  allowing  y  to  increase  after  the  limit  state.  The  general 
features  of  the  model  are  shown  in  Fig.  1,  which  also  illustrates  how  the 
presence  of  a  extends  the  elastic  regime. 

d 

An  appropriate  function  for  y  is  developed  by  assuming  that  a  path-length 
inelastic-strain  invariant,  e1 ,  exists  and  that  it  has  a  unique  value,  e^,  at 
the  limit  state  no  matter  which  path  is  followed.  One  function  that  provides 
the  initial  value  of  y  =  Yo  when  e*1  =  0  and  the  limit-state  value  of  y  =  yl 
when  e1*  e£  is  the  following: 


P 

Figure  1.  General  features  of  viscoplastic  model  in  L-P  plane. 


The  parameter  n  can  be  chosen  to  provide  theoretical  curves  that  have  the 
characteristic  elastic-plastic  shape  of  the  experimental  stress-strain  curves, 
although  detailed  fits  cannot  be  expected  when  merely  a  one-parameter  function 
is  used. 

Because  an  associated  flow  rule  provides  too  much  dilatation,  a  potential 
function  of  the  following  form  is  used: 

4>*  =  Y*P  -  L  -  a  (6) 


where 


Y*  «  y  +  b  ( 

and  b  is  a  material  constant.  A  nonassociated  flow  rule  implies  that  the 
inelastic  strain  increment  is  given  by 


i  36* 

de1  =  dx  g|- 


*  D  To 


C4-h) 


in  which  a  and  I  denote  the  stress  and  identity  tensors,  respectively,  and  dx 
is  a  positive  parameter  whose  value  is  obtained  by  solving  the  consistency 
condition 


$  s  0 


(9) 


Equation  8  shows  that  if  b  is  positive,  dilatation  is  reduced  from  that 
given  by  the  associated  flow  rule.  For  hydrostatic  compaction,  and  if  o$  is 
negligible  with  respect  to  P,  L  is  approximately  equal  to  P  and  the  volumetric 
inelastic  strain  increment  is  given  by 

de’  =  trace  dg1 

=  dX(l  -  y  -  b)  (10) 

The  requirement  that  a  hydrostatic  loading  path  should  always  yield  compaction 
is  satisfied  if  b  >  1  -  y^.  If  only  triaxial  loading  paths  are  available,  b 
is  chosen  by  trial  and  error  to  provide  an  optimal  fit  for  the  ensemble  of 
paths. 

For  metal  plasticity,  “e*  is  chosen  to  be  the  path  length  of  the  second 
invariant  of  the  inelastic  strain  deviator.  The  use  of  the  second  invariant 
for  soils  is  inappropriate  because  inelasticity  occurs  in  soils  subjected  to 
hydrostatic  loading.  It  has  been  found  (12)  that  a  weighted  combination  of 
first  and  third  invariants  is  useful  for  defining  a  strain  invariant  for  con¬ 
crete.  It  is  postulated  that  the  sane  approach  can  be  used  for  soils.  There¬ 
fore,  the  differential  of  the  path-length  strain  invariant  is  defined  by 

de1  (di!  +  cdi 31)  (11) 


in  which  c  and  P0  are  material  parameters,  and 

del  =  |  tr  de1  | 

dej  =  |  det  (dg1)  | 1/3 


(12) 

(13) 
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represent  differentials  of  the  absolute  values  for  the  first  and  third  invari¬ 
ants,  respectively.  If  P0  is  large,  the  effect  of  P  is  minimal.  If  P0  is 
small,  e1  accumulates  less  rapidly  as  P  increases,  a  feature  exhibited  by  many 
frictional  materials.  The  parameter  c  controls  the  relative  importance  of  the 
first  and  third  invariants. 

The  use  of  the  invariant  e1  is  useful  only  if  is  the  same  for  all 

possible  paths.  Inevitably,  there  will  be  scatter  in  the  experimental  value 
_i 

of  ei_  if  data  are  available  for  a  number  of  paths,  but  P0  and  c  can  be  chosen 
to  minimize  the  dispersion. 

If  linear  isotropic  elasticity  is  assumed,  the  specification  of  Young's 
modulus,  E,  and  Poisson's  ratio,  v,  completes  the  model  except  for  the  rate 
effects.  Because  very  few  data  involving  strain  rate  are  available,  the  fol¬ 
lowing  simple  relation  is  assumed: 


The  total  strain  invariant,  e,  is  defined  by  using  the  same  parameters  as 
those  used  for  "C1  in  Eq.  11.  The  parameters  ao  and  m  are  chosen  to  provide 
the  general  characteristics  of  enhancement  as  indicated  by  experimental  data. 

The  reference  strain  rate,  ¥0,  is  chosen  merely  to  make  the  argument  dimen¬ 
sionless  for  the  exponential  operation. 

The  model  outlined  in  this  section  follows  the  concepts  of  classical 
plasticity;  therefore,  the  development  of  an  algorithm  for  incorporating  the 
constitutive  equation  into  computer  codes  presents  no  problem.  The  material 


parameters  E,  v,  o$,  Y|_,  n,  b,  c,  P0,  o0,  and  m  represent  a  large  but  manage¬ 
able  set,  especially  In  light  of  the  wide  range  of  response  characteristics 
that  can  be  predicted.  In  the  next  section  material  parameters  and  theoreti¬ 
cal  predictions  are  given  for  the  response  of  two  sands  under  both  static  and 
dynamic  conditions. 

COMPARISONS  WITH  EXPERIMENTAL  DATA 

Large  strain-rate  testing  of  soils  presents  many  problems  for  the  exper¬ 
imentalist.  For  example,  if  a  specimen  is  loaded  too  rapidly,  the  propagation 
of  stress  waves  will  create  nonuniform  conditions  of  stress  and  strain  in  the 
sample.  Whitman  (13)  points  out  that  In  the  case  of  triaxial  samples,  the 
effect  of  lateral  Inertia  must  be  Included  in  a  comprehensive  assessment  of 
failure  and  material  properties.  For  specimens  4  inches  long  and  1.5  inches 
in  diameter,  he  notes  that  the  time  to  failure  should  be  greater  than  5  ms  if 
lateral  inertia  is  to  be  ignored. 

In  the  case  of  dynamic  uniaxial  strain  testing,  the  rate  of  loading  is  so 
rapid  that  no  drainage  occurs,  particularly  In  fine-grained  soils.  Here  thin¬ 
ner  samples  are  typically  used  in  order  to  prevent  the  formation  of  nonuniform 
conditions  of  stress  and  strain  as  well  as  to  reduce  side  friction  effects 
between  the  sample  and  equipment  boundaries.  However,  Schindler  (10)  has 
developed  equipment  that  can  be  used  to  test  samples  as  large  as  10  inches  in 
diameter  and  2-1/2  inches  thick.  The  advantages  of  testing  larger  samples  are 
reported  to  be  (1)  reduced  disturbance  at  the  boundaries,  (2)  increased  magni¬ 
tude  of  the  displacements  recorded  and  thereby  improved  accuracy  of  the  mea¬ 
surements,  and  (3)  measuranent  of  displacement  at  several  points  on  the  sample 
face.  It  can  be  shown  that  if  the  rise-time  of  the  loading  pulse  is  long 


relative  to  the  time  required  for  the  propagating  stress  wave  to  traverse  the 
distance  back  and  forth  between  a  rigid  boundary  and  a  free-surface  boundary, 
then  inertial  stresses  may  be  neglected.  This  principle,  called  the  multiple 
reflection  technique,  has  been  the  basis  for  the  construction  of  many  dynamic 
uniaxial  strain  test  facilities. 

The  triaxial  test  apparatus  used  by  Ito  and  Fujimoto  (3),  in  which  an 
hydraulic  oil  pump  was  employed,  was  reported  to  be  capable  of  providing  axial 
loads  with  variable  displacement  rates  of  from  0.04  mm/s  to  600  mm/s.  A  load 
cell  was  used  to  measure  axial  loads  and  an  eddy  current  type  transducer  to 
measure  axial  displacements.  Radial  displacements  were  recorded  by  six  opti¬ 
cal  displacement  transducers  having  a  sensitivity  of  0.01  mm.  The  average 
radial  displacement  was  used  in  the  volume  change  calculations.  Air  was  used 
to  apply  the  confining  pressure.  In  order  to  reduce  friction  at  the  sample 
ends,  two  sheets  of  rubber  lubricated  with  silicon  grease  were  placed  at  each 
end  of  the  sample,  between  the  sample  and  the  end-plate. 

Two  types  of  uniaxial  strain  test  equipment  were  used  by  WES  (4).  The 
first  was  a  100-MPa  gas-driven  ram  device,  similar  in  concept  to  the  apparatus 
described  by  Schindler  (10).  Capabilities  include  piston  loads  up  to 
450,000  N,  rise-times  as  short  as  3  ms,  and  decay  times  as  short  as  20  ms. 

The  piston  pressurizes  an  oil-filled  chamber  that  produces  the  load  on  the 
sample  face.  Because  of  the  high  lateral  stresses  created  in  the  sample,  the 
device  must  be  very  rigid  in  order  to  prevent  radial  deformations  from  occur¬ 
ring.  In  the  second  device,  recently  developed  by  WES,  an  explosive  charge 
drives  the  piston  into  the  oil.  Pressures  up  to  100  MPa  may  be  obtained,  with 
rise-times  and  decay  times  as  short  as  0.3  ms.  The  samples  used  in  these 
devices  are  quite  thin  (1.27  cm)  because  of  the  previously  mentioned  require¬ 
ment  that  the  rise-time  of  the  loading  be  long  in  relation  to  the  transit  time 


if  the  multiple  reflection  theory  is  to  be  valid.  The  sample  height-to- 
diameter  ratio  was  1:76  in  the  WES  tests. 

For  the  data  provided  by  Ito  and  Fujimoto  (3),  the  model  variables  were 
adjusted  to  obtain  an  adequate  fit  to  the  static  triaxial  compression  tests. 
The  variables  used  were  E  =  1300  MPa,  v  =  0.25,  e^  =  0.19,  <j$  =  0.007  MPa, 

Y0  =  1.06,  y ^  =  0.86,  n  =  0.45,  b  =  0.33,  c  =  1.8,  P0  =  -5  MPa,  a0  =  0.19  MPa, 
and  m  =  0.75.  Two  tests,  conducted  at  different  confining  pressures  (0.2  and 
0.3  MPa),  were  available  for  defining  parameters.  The  samples  of  Toyoura  sand 
used  in  these  tests  were  air-dried  and  had  an  initial  void  ratio  of  0.72.  The 
experimental  results  and  theoretical  predictions  for  the  static  tests  are 
shown  in  Fig.  2  as  plots  of  octahedral  shear  stress  (T0)  versus  octahedral 
shear  strain  (c0).  In  terms  of  principal  stresses  and  strains,  t0  and  c0  are 
given  by 

T0  =  -J  [«L  -  «2)2  ♦  -  O3)2  ♦  ("3  -  =,)2]1/2  (15) 

Co  »  -J  C(«1  -  e2)2  +  <e2  -  e,)2  ♦  (e3  -  ei)2]1/2  (16) 

Corresponding  results  are  shown  in  Fig.  3  for  the  case  in  which  the  axial 
strain  rate  is  1/s.  Both  the  theoretical  and  the  experimental  results  show 
some  enhancement  of  strength  together  with  a  shift  in  strain  at  the  limit 
state.  The  model  predicts  a  steeper  response  up  to  a  higher  level  of  stress 
for  the  rate-effect  case  than  for  the  static  case.  This  prediction  is  sup¬ 
ported  by  experimental  data  when  the  confining  pressure  is  0.3  MPa  but  not 
when  the  confining  pressure  is  0.2  MPa.  The  reason  for  this  discrepancy  is 
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(b)  Confining  pressure  0.3  MPa 

Figure  2.  Theoretical  and  experimental  (3)  results 
for  triaxial  compression  of  sand. 
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(a)  Confining  pressure  0.2  MPa. 


(b)  Confining  pressure  0.3  MPa. 

Figure  3.  Theoretical  and  experimental  (3)  results 

for  triaxial  compression  of  sand,  longitudinal 
strain  rate  1/s. 


Volumetric  strain  as  a  function  of  octahedral  strain  for  the  triaxial 
path  with  a  confining  pressure  of  0.3  MPa  is  shown  for  the  static  and  dynamic 
cases  in  Figs.  4a  and  4b,  respectively.  No  volumetric  data  for  a  confining 
pressure  of  0.2  MPa  were  presented  in  Ref.  3.  Slightly  more  dilatation  is 
predicted  for  the  dynamic  case  than  for  the  static  case,  but  the  experimental 
data  show  a  significant  reduction  in  dilatation.  Whitman  (13)  indicates  that 
more  rather  than  less  dilatation  should  be  expected  with  higher  rates  of  load¬ 
ing,  but  there  is  no  experimental  evidence  to  corroborate  this  statement. 

The  stress-strain  results  illustrated  in  Figs.  2  through  4  represent  the 
response  of  sand  subjected  to  relatively  small  confining  pressures  and  strain 
rates.  The  WES  data  (4),  on  the  other  hand,  represent  the  response  of  a 
clayey  sand  subjected  to  confining  pressures  higher  by  an  order  of  magnitude 
and  strain  rates  higher  by  two  orders  of  magnitude  than  those  imposed  upon  the 
sand.  The  representation  of  both  loading  situations  by  one  model  indicates 
that  fundamental  response  features  are  contained  within  the  constitutive 
relation. 

The  following  parameters,  chosen  to  fit  the  experimental  data  for  the 
clayey  sand,  were  used  in  the  analysis  of  the  WES  data:  E  =  1200  MPa, 
v  =  0.25,  o$  =  0.1  MPa,  e^  =  0.55,  yo  a  1.2,  Y|_  =  0.82,  n  =  0.4,  b  =  0.35, 
c  =  3.2,  P0  =  -32.5  MPa,  o0  =  0.007  MPa,  and  m  =  0.30.  The  data  consisted  of 
stress  and  strain  results  for  two  static  triaxial  conpression  tests  performed 
at  confining  pressures  of  2.07  MPa  and  4.14  MPa  and  for  a  uniaxial  strain  test 
during  which  lateral  stress  measurements  were  made.  Figs.  5  and  6  compare  the 
theoretical  predictions  with  the  experimental  data.  The  agreement  is  consid¬ 
ered  excellent. 
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.  Theoretical  and  experimental  (3)  volumetric  results 
for  triaxial  compression,  confining  pressure  0.3  MPa. 
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The  WES  dynamic  data  were  obtained  from  dynamic  uniaxial  strain  tests, 
each  of  which  was  performed  at  varying  strain  rates.  The  strain  histories  for 
two  cases  (Fig.  7)  were  obtained  by  choosing  strain  and  time  values  from  plots 
of  vertical  stress  versus  time  and  vertical  stress  versus  vertical  strain 
presented  in  Ref.  4.  The  maximum  uniaxial  strain  rate  is  approximately  200/s 
for  both  cases. 

Spline  fits  to  the  experimental  data  (Fig.  7)  were  used  as  input  to  the 
model  to  obtain  theoretical  values  of  stress.  Fig.  8  is  a  comparison  of  the 
theoretical  predictions  and  experimental  results  for  the  dynamic  uniaxial 
tests.  The  response  of  the  soil  is  extremely  different  from  that  obtained  in 
the  static  test.  Here  the  vertical  stress-vertical  strain  curve  is  convex  to 
the  strain  axis,  whereas  the  static  test  shows  the  more  conventional  response: 
concave  toward  the  strain  axis.  The  viscoplastic  model  has  captured  this 
feature  at  least  qualitatively. 

The  softening  response  is  due  to  a  decrease  in  strain  rate  rather  than  to 
the  actual  strain-softening  type  of  behavior  that  might  be  seen  in  static 
triaxial  tests.  Unfortunately,  lateral  stress  measurements  were  not  measured 
for  the  dynamic  test,  and  no  dynamic  triaxial  test  results  are  available  for 
this  material . 

CONCLUSIONS 

The  elasto-viscoplastic  model  presented  here  has  been  shown  to  be  capable 
of  predicting  the  stress-strain  and  volumetric  strain  response  of  soils  (sand 
and  clayey  sand)  to  a  selected  group  of  tests  (triaxial  compression  and  uniax¬ 
ial  strain).  The  model  has  been  used  to  predict  soil  response  to  strain  rates 
as  high  as  200/s.  The  dramatic  changes  in  soil  behavior  associated  with  rate 
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Vertical  strain  as  a  function  of  time  for  WES  data 
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effects  are  reproduced  by  the  model.  The  combination  of  mathematical  simplic¬ 
ity  and  realistic  predictions  make  the  model  particularly  appropriate  for 
computer  codes  used  to  study  soil  response  to  dynamic  loading. 
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NOTATION 


The  following  synbols  are  used  in  this  paper: 

b  *  material  parameter; 
c2  3  material  parameter; 
dx  =  path-length  parameter; 

E  =  Young's  modulus; 
ei»  e2»  e3  =  principal  strains; 

e^  =  inelastic  volumetric  strain; 

J  3  total  strain  invariant; 
e^  =  path-length  inelastic-strain  invariant; 

*  path-length  strain-invariant  limit  value; 

V0  *  reference  strain  rate; 
g1  3  inelastic  strain  tensor; 


l  3  identity  tensor; 

L  3  stress  invariant; 


m  =  material  parameter; 
n  *  material  parameter; 

P  «  mean  pressure; 

P0  *  material  parameter  with  units  of  pressure; 

Y  =  slope  of  flow  surface  in  L-P  space; 

Y ^  =  slope  of  limit  surface  in  L-P  space; 

Y0  ■  initial  slope  of  flow  surface  in  L-P  space; 
Y*  =  slope  of  potential  surface  in  L-P  space; 

50  *  octahedral  shear  strain; 
v  =  Poisson's  ratio; 

a  -  additional  cohesion  due  to  strain  rate; 
d 

as  =  cohesion; 
a0  =  material  parameter; 
lt  o2,  o3  *  principal  stresses; 
a  »  stress  tensor; 

*0  =  octahedral  shear  stress; 

4  =  flow  function;  and 
*  potential  function. 
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